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Preambule 



L'objectif de ce cours est de presenter et d'appliquer (TD) les outils de base qui permettent 
de decrire et predire les phenomenes de transfert de matiere, de quantite de mouvement (forces) et 
d'energie. De tels phenomenes ont lieu non seulement dans des equipements construits par 
rhomme mais egalement dans le milieu naturel, de l'echelle de la cellule jusqu'a celle de la 
biosphere. 

La matiere qui nous entoure couramment possede une tres grande densite moleculaire 
(nombre de molecules par metre cube) on peut sans probleme la decrire a partir des concepts 
introduits en thermodynamique et la considerer comme un milieu continu. De ce fait, on peut tirer 
profit d'un grand nombre d'outils de l'analyse et de la geometric 

Les transferts de matiere, de quantite de mouvement et d'energie presentent des similitudes, 
on presentera un formalisme unique pour ecrire les differents bilans que ce soit sur des systemes 
macroscopiques (equations differentielles ordinaires) ou des volumes elementaires (equations aux 
derivees partielles). 

Les travaux diriges permettront d'aborder un certain nombre d'applications de ces outils et 
les etudiants auront largement l'occasion d'utiliser et de completer ces bases dans des 
enseignements optionnels. 

Ce cours s'inspire essentiellement de l'ouvrage " Transport Phenomena " de R.B. Bird, 
W.E. Stewart et E.N. Lightfoot (John Wiley & Sons) ou les etudiants trouveront de nombreux 
complements et diverses illustrations. On pourra egalement se reporter a l'ouvrage "Mecanique et 
Rheologie des fluides en Genie chimique" de N.Midoux (Tec. Doc. Lavoisier). 
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Introduction 



Une seconde methode d' application de la mecanique, qui evite de prendre en consideration la division de la matiere 
jusqu 'aux points materiels « reels », est la mecanique des soi-disant milieux continus. Cette mecanique est 
caracterisee par la fiction que la densite et la vitesse de la matiere dependent d'une fagon continue des coordonnees 
et du temps, et que la partie des interactions non explicitement donnee peut s' interpreter comme des forces 
superficielles qui sont aussi des fonctions continues de la position. Albert Einstein 

La physique est toujours un petit peu fausse Yves Rocard 

La mecanique des milieux continus considere que la matiere remplit continument des domaines 
determines de l'espace eventuellement separes par des interfaces. Ainsi a la pression 

17 

atmospherique, le volume moyen occupe par une molecule de gaz parfait est d'environ 4.10" 
mm . II est possible, sur des volumes tres petits devant l'echelle macroscopique, de definir une 
valeur moyenne d'une propriete de ces particules qui soit intrinseque, au sens par exemple, ou la 
moyenne effectuee sur une moitie des particules est egale (a un infiniment petit pres) a la moyenne 
effectuee sur l'autre partie. L'exemple le plus connu d'une de ces proprietes moyennes 
macro scopiques est l'energie interne d'un gaz monoatomique: il s'agit de l'energie cinetique des 
particules (lorsque celui-ci est au repos en moyenne macroscopique). 

L'objet de ce cours est de presenter et d'utiliser les outils qui permettent de decrire et de predire 
revolution (dans l'espace et dans le temps) de cette matiere a l'echelle macroscopique. Nous 
ferons largement appel a la thermodynamique pour: 

determiner les relations qui lient les differents parametres (p,p,T, H ...) d'un materiau 
homogene: equations d'etat. On supposera (hypothese de l'equilibre local) que dans un tres 
petit volume autour d'un point, la matiere peut etre consideree comme homogene. 

exemple : gaz parfait :p/p=RT/M, H{T} = h{t o } + C p . (T - T ) 

caracteriser l'equilibre lorsqu'il y a coexistence de deux etats possibles de ce materiau 
(equilibre liquide-vapeur, equilibre chimique) 

exemple : equilibre liquide-vapeur pour un corps pur : p=psat(T) 

formaliser les principes fondamentaux de conservation de la masse et de l'energie, 
auxquels nous ajouterons la relation fondamentale de la mecanique, c'est-a-dire le principe 
de conservation de la quantite de mouvement. 

accumulation + bilan des flux (sortie-entree) = terme source 

quantifier les flux diffusifs de matiere et de chaleur, auxquels il faut rajouter ceux de 
quantite de mouvement qui correspondent aux frottement visqueux. Nous nous limiterons a 
l'approximation lineaire (non couplee) qu'offre la thermodynamique des phenomenes 
irreversibles. 

exemple : flux diffusif de chaleur oc gradient de temperature 
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1) Transfert de matiere 



Nous allons commencer par analyser le transport d'un composant (dilue, du sel par exemple) dans 
un fluide (de l'eau par exemple) dont le mouvement est connu. Dans un premier temps nous 
considerons un petit element de volume parallelepipedique (Ax.Ay.Az) fixe. 

1.1) Equation de conservation sur un element de volume fixe (masse totale et constituants) 
1.1.1) Masse totale 

Nous savons que la masse totale se conserve. Du fluide rentre ou sort par chacune des 6 facettes de 
l'element de volume en fonction de la vitesse. 



v{x + Ax/2,y,z} 




^ v x {x + Ax/2,y,z} 



x-Ax/2 



x+Ax/2 



L'accumulation de masse dans l'element entre t- At/2 et t+At/2 vaut: 

[ p{x,y,z,t+At/2} - p{x,y,z,t-At/2} ]. Ax.Ay.Az (developpement limite a l'ordre 2 en x,y,z et t) 

La masse qui sort par la facette situe a Vest entre t-At/2 et t+At/2 vaut: 
p { x+Ax/2,y,z } . v x { x+Ax/2,y,z } . Ay.Az. At 

Ce flux est du au mouvement du fluide, on dit qu'il s'agit d un flux convectif. 

Si on divise le bilan sur l'element par le volume de l'element on obtient alors: 



accumulation 



p{ x, y, z, t + At / 2} - p{x, y, z, t - At / 2} 



At 



(p. v x ){x + Ax / 2, y, z, t} - (p. v x ){x - Ax / 2, y, z, t} 



bilan des flux 
convectifs 



Ax 



(p.v y ){x,y + Ay/2, z, t}-(p.v y ){x,y- Ay/2,z,t} 
Ay" 

(p.v z ){x, y,z + Az/2, t}-(p.v z ){x,y,z-Az/2, t} 
Az 



= 
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Si Ton fait tendre Ax, Ay, Az et At vers zero (les termes d'ordre plus eleve des developpements 
limites disparaissent), on obtient: 



3p 3(p.v x ) 3(p.v y ) 3(p.v z ) 

3t 3x 3y 3z 



accumulation bilan (sortie - entree) des flux convectifs 



Si on note pv 



vpv z ; 



le vecteur densite de flux convectif 



3x 
_3_ 

3y 
_3_ 

^3z^ 



le vecteur symbolique nabla (ici en coord, cartesiennes) 



alors cette relation peut encore s'ecrire: 



+ V.(pv) = 
accumulation bilan des flux convectifs 



L'expression V»( f\) encore notee Div(fj)est appelee divergence du vecteur densite de flux f\. 
II s'agit du bilan (sortie moins entree) du flux qui traverse les facettes de l'element, divise par le 
volume de cet element. Remarquez la presence du point qui represente le produit scalaire de deux 
vecteurs. 



1.1.2) Masse du constituant i 



La meme demarche peut etre utilisee pour le bilan de la masse d'un constituant i sur un element de 
volume. Notons 0)j la fraction massique du constituant i (masse du constituant i / masse totale). La 
masse du constituant i par unite de volume est alors p(0[ . Laccumulation de masse du constituant i 
par unite de temps et de volume vaut done: 

3pco ; 
~3T 

Le flux du constituant i a travers les facettes de l'element peut etre du au mouvement du fluide 
qui transporte avec lui ce constituant. Ce flux est appele convectif. 

Le vecteur densite de flux convectif du constituant i vaut : p(0 i v 

Sa direction et son sens sont celui du vecteur vitesse. Sa norme s'exprime en (kg{i}/s) /m 2 . 
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Le flux du constituant i peut egalement etre du au fait que la fraction massique Q0[ n'est pas 
homogene dans le milieu. Ce flux est appele diffusif. 

Le vecteur densite de flux diffusif du constituant i est note: j; . 

II depend du gradient de la fraction massique (ou plus generalement du gradient du potentiel 

electrochimique). Son sens est du plus concentre au moins concentre. Sa norme s'exprime 
egalement en (kg{i}/s) / m 2 . 




Le vecteur densite de flux total du constituant i s'ecrit simplement: fj; = pa>jV + j ; 

Le bilan des flux du constituant i a travers les facettes de l'element de volume divise par le volume 
de l'element s'ecrit : V.f); . 

Finalement l'equation de conservation de la masse du constituant i s'ecrit: 



+ V. (p©jV + jj) = r, 

accumulation bilan des flux convectifs terme source 

et diffusifs 



Le terme r; est la production du constituant i par unite de temps et de volume (kg{i}/s/m ) du fait 
d'eventuelles reactions chimiques. 

Remarque: 

La somme des densites de flux de tous les constituants vaut par definition pv , d'ou: 
pv = S f li =P(£ ( ° i )v + Zl orEcOi=l ^Xl =0 

i 1 i 1 

C'est a dire que la somme des flux diffusifs de tous les constituants est nulle. 

Si Ton somme les equations de conservation de la masse de tous les constituants on retrouve 
l'equation de conservation de la masse totale. 

La vitesse v est celle du barycentre de tous les constituants, c'est-a-dire la vitesse moyenne des 
molecules ponderee par leurs masses. 
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1.2) Equation de conservation sur une particule de fluide suivie dans son mouvement . 
(derivee particulaire) 

Les equations de conservation de la masse totale et de celle d'un constituant i s'ecrivent: 

3p 3(p.vJ ^ 3(p.v y ) 3(p.v z ) = 

3t 3x 3y 3z 

3pa); 3(pQ) i .v x ) 3(pQ>i.v y ) 3(p(Dj.v z ) 3j ix 3j iy 3j iz 

IT + ^x~ + + ~^z— + "aT + 17 + "3z" = r; 



en multipliant la premiere equation par co ; et en la soustrayant a la seconde on obtient: 
(enutilisant(f.g)'=f.g + f.g') 



3(0; 

'~3T 



+ pv x 



3(0; 

3x~ 



+ p V 



3(0; 

17 



+ p v 2 



3(0; 

~3z~ 



3j ix , a J.y 3j iz 
+ — + — + — = r 



3x 3y 



3z 




Le terme entre crochet est appele derivee particulaire de (Oj , on le note 



D(0; 

~Dt 



La derivee particulaire de (0; represente revolution (au cours du temps) de 0)i lorsque Ton suit une 
particule de fluide (de masse totale fixee) dans son mouvement. 

Sous forme vectorielle l'equation devolution de la fraction massique en suivant le fluide dans son 
mouvement s'ecrit : 



D(0; 

Dt 



3(0; 



dt 



L + (V.V)(0; 



Remarque : la particule de fluide suit le mouvement et se deforme avec lui, il n'y a done pas de 
mouvement moyen a travers ses frontieres et done pas d'echange convectif. Elle n'echange de la 
matiere avec l'exterieur que par diffusion. 

1.3) Equation d'etat 

Dans le cas general, la masse volumique est une fonction de la fraction massique des differents 
constituants (ainsi que de la pression et de la temperature) : p = f {p, T, (0; } 

On a cependant souvent affaire a des milieux dilues, dans ce cas on peut considerer que la masse 
volumique est constante. 

Remarquez que le volume ne se conserve pas a priori: un litre d'eau plus un litre d'alcool ne 
forment pas deux litres de melange. 
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1.4) Equation de transfert (loi de Fick) 



1.4.1) Transfert unidirectionnel 

Considerons un melange binaire non homogene de molecules A • et B O au repos d'un point de 
vue macro scopique. Supposons que la fraction massique en A n'est fonction que de y et diminue 
avec y , ce qui se traduit par : dO) A /dy < 0. 

L'agitation thermique ( O ^ ) fait qu'une molecule a une chance sur deux d'aller vers les y 
positifs ou negatifs. Si Ton fait le bilan (a travers un plan perpendiculaire a oy), on trouve un flux 
net du constituant A dans le sens des y positifs, ce qui se traduit par : j A . y >0 . 

AA A 
D O • o o 

A A A 2 particules • montent et une seule descend 



V V v 

» • o o • • 

V v V 

d(0 A 

Un modele simple qui puisse traduire cela est celui de Fick: j A = -(constante) — — 

La loi de Fick est tres largement utilisee, en phase gazeuse ou liquide, pour estimer les flux 
diffusifs. Elle peut etre generalisee comme nous l'avons vu en thermodynamique en utilisant la 
variation du potentiel electrochimique au lieu de celui de la fraction massique. La diffusivite est en 
general fonction de la pression, de la concentration et de la temperature (elle augmente avec la 
temperature) 

d(0 A r i 

j A =- D A B p-^ avec D AB = f|(Q A ,p,Tj 

1.4.2) Loi de transfert pour un champ quelconque de concentration 

Nous avons suppose precedemment que la fraction massique etait lineaire par rapport a y (et 
independant de x). La surface (0^=f{x,y}est dans ce cas particulier un plan dont la pente est nulle 
dans la direction ox et dont la plus grande pente est dirige selon oy. Le vecteur densite de 
flux diffusif est oriente selon -oy, c'est a dire dans la direction de la plus grande pente 
descendante, et la densite de flux est proportionnelle a cette pente. 
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Considerons maintenant que la fraction massique est une fonction quelconque de x et de y. La 
surface 0) J \=f{x,y} a maintenant un relief quelconque. En un point (x,y) nous pouvons trouver le 
plan tangent a cette surface, le vecteur densite de flux suit comme precedemment la ligne de plus 
grande pente descendante et est proportionnel a cette pente. 



A 



CO. 




Le vecteur oriente selon la plus grande pente ascendante et dont la norme est egale a cette pente est 
le vecteur gradient, note: 



Vo) A = Grad (0 A 



dx 

3® A 

V dz J 



Noter que cette fois le vecteur nabla V n'est pas suivi d'un point, il ne s'agit pas d'un produit 
scalaire, C0 A est un scalaire et V co A est un vecteur. 



Finalement, pour un melange binaire, la loi de Fick ,qui exprime le vecteur densite de flux 
diffusif s'ecrit: 



j A =-PDab V 0)a o 



f ■ \ 

Ja.x 



)A.y 



VJa.z^ 



= - P D 



AB 



3o) 



dx 
dy 

^A 

V dz ) 



Rappelons que par definition (0^ +(Oq =1 , Dab= Dba , j A = — j B 
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1.5) Conditions aux limites 



Les equations que nous venons de presenter permettent de decrire ce qui se passe dans un 
domaine rempli d'un milieu continu. Pour predire les evolutions du systeme il est indispensable de 
connaitre egalement ce qui se passe sur les frontieres du domaine, encore faut-il definir tres 
clairement les limites de ce domaine. 

Les conditions aux limites les plus simples peuvent se classer en trois categories qui portent 
chacune un nom celebre : 

Dirichlet : on connait les fractions massiques des constituants 
Neumann : on connait les flux diffusifs des constituants 

Cauchy : on connait une relation entre les fractions massiques et les flux des constituants 

a) Condition de Dirichlet : on connait la fraction massique 1 

C'est typiquement le cas a l'entree d'un systeme. On connait la composition a l'entree d'un 
reacteur, la charge en polluant a l'entree d'un digesteur, la teneur en sel d'une riviere alimentant 
un marais,etc. 

C'est souvent le cas s'il y a des equilibres de phase. A l'interface entre des cristaux de sel et une 
solution, la concentration est egale a la concentration de saturation (loi d'equilibre entre la forme 
cristalline et la forme solutee). Si on s'interesse a la diffusion de la vapeur d'eau au-dessus d'un 
etang, on peut dire qu'a l'interface la pression partielle de vapeur d'eau est egale a la pression de 
vapeur saturante (qui est fonction de la temperature de l'eau a l'interface : loi d'equilibre liquide- 
vapeur pour un corps pur) 

Cela peut etre une bonne approximation a l'interface avec un milieu agite. Considerons par 
exemple un milieu poreux sature en eau qui est plonge dans une solution aqueuse agitee et de 
concentration connue. Au niveau de l'interface, la concentration de la solution contenue dans le 
milieu poreux est proche de celle du liquide. Dans ce cas, la diffusivite du solute a travers l'eau 
contenue dans le milieu poreux a une signification et une valeur differente de la diffusivite du 
couple solute/eau 2 . 

b) Condition de Neumann, on connait la densite de flux diffusif du constituant 

C'est typiquement le cas pour une paroi impermeable. La densite de flux total (convectif + 
diffusif) est nulle, par ailleurs, la vitesse est nulle pres d'une paroi il n'y a done pas de flux 
convectif. On en deduit que la densite de flux diffusif est nulle. 

On peut souvent considerer qu'a la sortie d'une dispositif la diffusion est negligeable devant la 
convection 3 . 



Si le milieu est un melange de I constituants, il faut imposer les fractions massiques pour (1-1) constituants, 
la fraction massique du Ieme constituant etant obtenue par la relation E (D i = 1 

2 Si la matrice du milieu poreux interagit avec le constituant on peut plus ecrire Feaglite des concentration mais 
seulement l'egalite des potentiels electrochimiques du constituant des deux cotes de l'interface 

3 II faut pour cela que la vitesse en sortie soit suffisamment importante, la tuyauterie assez longue jusqu'au 
prochain dispositif, ceci compare a la diffusion moleculaire ou turbulente, v.L/D 
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D'autres situations peuvent etre approchees par ce type de condition. Si, par exemple, on envoie 
une solution tres concentree avec un debit tres faible a travers une paroi poreuse dans un systeme 
ou la solution est peu concentree. On connait done la densite de flux total du constituant. Du cote 
de la solution, la vitesse est tres faible, le transfert s'y fait done essentiellement par diffusion. Tout 
se passe done comme si on connaissait la densite de flux diffusif a l'interface du cote de la 
solution. 



vitesse quasi nulle — > transfert par diffusion 

(sous l'effet du gradient de concentration) 




densite de flux milieu domaine d'etude 

connue poreux des transferts de matiere 



Dans ces differents cas, on connait done la densite de flux diffusif , que nous noterons ji. en trant , a 
travers une section ; comment traduire cela en equation ? 



Commencons par une section 
perpendiculaire a ox 



y 



Dans ce cas nous connaissons en fait la composante selon ox de la densite de flux diffusif . 
que nous pouvons exprimer par ailleurs par la loi de Fick : 



Ji.entrant 



3(0; 

ax" 



Si la section est perpendiculaire a oy , on a 

aco : 



Ji, 



= j'i. y =-pDi 



a y 



Dans le cas d'une section perpendiculaire a un vecteur n (vecteur normal unitaire) on connait en 
fait la projection sur n du vecteur densite de flux diffusif. On obtient cette projection en prenant 
le produit scalaire du vecteur densite de flux diffusif j ; et de n . 



Comme, par convention n est oriente vers l'exterieur du domaine, ce produit scalaire est egal au 
flux qui sort, e'est-a-dire a l'oppose de ji.entrant • Finalement : 
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Exemple : paroi inclinee a 45 degres 
n = (-1/ 4l -1/ 4l) 



interieur du domaine 



Ji.e 



pD^Vco^.n 



= -pD i 



1 3(0, 1 3<o ; 
+ 



V2 3x V2 By 




Si la paroi est impermeable on aura done : j ientrant = pD ; (V(0j).n = => (V(O ; ).n = 

1 3(0; 1 3(0; 

pour une paroi a 45 degres : — + — = ^ 

On verifiera que pour des sections perpendiculaires a ox ou oy , on retrouve bien les relations 
precedentes. 

c) Condition de Cauchy : on connait une relation (affine) entre la densite de flux 

diffusif du constituant et la fraction massique 

Donnons un exemple. La frontiere est constitute d'une membrane de faible epaisseur: e (qui 
n'accumule pas de matiere), du cote exterieur on a une solution de concentration connue et 
fortement agitee. Sur la surface exterieure de la membrane, la concentration est alors proche de la 
concentration moyenne du milieu agite : (Ouxt • 



interieur du domaine 



(Oi.ext 



milieu exterieur 
parfaitement agite 



concentration 
uniforme 



OOOOO 
OOOOO 
OOOOO 
OOOOO 



(Oi.s 



gradient de la fraction massique 
perpendiculairement a la surface 



Si la vitesse moyenne est nulle dans la membrane alors les flux sont uniquement diffusifs (pas de 
convection puisque v = dans la membrane et a 1' interface entre la membrane et le domaine 
etudie) 

Le flux traversant la membrane (par diffusion) s'ecrit 

C0 iex( -(0; S 

ji entrant = P D * — " '~ = P- 1 (C0;. ext - (O i s ) 



Dans cette expression, D* est la diffusivite du solute a travers la membrane (cette diffusivite a une 
signification et une valeur differentes de la diffusivite du constituant i dans de l'eau seule). Le 
coefficient k est appele coefficient de transfert de matiere de la membrane. 
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Par ailleurs, le flux qui traverse la membrane est egal au flux diffusif a l'interface qui s'ecrit : 

ji.entrant = P D i V©; | g . fi 

A partir de ces deux expressions on obtient alors une relation entre la fraction massique et le 
gradient de la fraction massique a la limite du domaine 4 : 

P-k(co, ext -co iS ) = pD i Vco i | s .n 



pour une surface perpendiculaire a ox, on a 




Certaines membranes peuvent egalement avoir des affinites (sur l'une des faces) pour des ions de 
signe donne, par ailleurs la vitesse moyenne n'est en general pas nulle a travers la membrane. La 
situation devient alors plus complexe a modeliser. 

D'une maniere generale, de nombreux modeles supposent que le flux a travers une membrane est 
proportionnel a la difference de fraction massique de part et d' autre de la membrane. Attention : 
parfois on se base plutot sur des differences de concentration ou de pression partielle dans ce cas le 
coefficient de transfert s'exprime dans des unites et prend des valeurs differentes (bien verifier 
dans quel cas on se trouve). 

Ji.entrant = P- k -(<0, e xt 0U ji.entrant = k *(C,ext "C^) OU j, e „ = k * * (p.^ - p.^ ) 



Generalites sur les conditions aux limites 

Dans la pratique les conditions aux limites ne peuvent pas toujours etre mises sous l'une des trois 
formes a-b-c precedentes, on peut cependant enoncer la regie suivante: Si Ton traite d'un melange 
de I constituants, il faut ecrire (1-1) equations scalaires portant sur les fractions massiques, sur 
l'ensemble de la frontiere. La relation suivante relie en effet les fractions massiques entre elles: 
E tOj = 1 

D'un point de vue mathematique ces equations doivent relier les variables intensives (fractions 
massiques) et les gradients (perpendiculairement a la surface) de ces variables intensives. 

D'un point de vue physique elles relient les variables transportees (fractions massiques) et les flux 
diffusifs et convectifs des quantites transportees. 

Nous verrons des conditions aux limites tout a fait similaires en transfert de chaleur, la variable 
intensive etant alors la temperature. 

Recapitulatif : 

Dirichlet : fraction massique imposee : (C0i)si connu 
Neumann : densite de flux diffusif imposee : (vo);)^.!! connu 

Cauchy : coefficient de transfert k donne : (vo); ) .ii = a.(a); ) sj + b , a et b connus 



Dans ce cas, a la limite du domaine, le flux diffusif pDjVcOj.n est une fonction affine de la fraction massique tOj s 
fraction massique a l'exterieur co i ext etant une donnee du probleme). 
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Conditions d' interface 



Dans certains des exemples de conditions aux limites vus precedemment, nous avons suppose que 
les caracteristiques du milieu exterieur au domaine d'etude etaient connues, y compris au niveau 
des interfaces. Si ce milieu exterieur est parfaitement agite, les variables intensives telles que la 
fraction massique peuvent etre considerees comme uniformes, leurs valeurs a l'interface sont alors 
proches des valeurs moyennes qui peuvent etre plus facilement connues. 

Parfois il n'est pas possible de faire de telles hypotheses, il faut alors traiter un systeme de deux 
milieux continus separes par une interface. Cette interface peut etre la surface de separation de 
deux phases (liquide et vapeur, par exemple) ou une paroi qui laisse passer un flux de la quantite 
transportee (membrane pour le transfert de matiere, paroi metallique d'un echangeur pour le 
transfert de chaleur). 



1.6) Combinaison des equations - Adimensionnalisation 
1.6.1) Combinaison des equations 

Nous allons maintenant combiner les differentes equations dans le cas particulier d'un melange 
binaire A/B, de masse volumique constante, repondant a la loi de Fick, sans reaction chimique, 
avec des conditions aux limites donnees par un coefficient de transfert uniforme. 

Prenons l'exemple schematise ci-dessous et raisonnons en deux dimensions. II pourrait s'agir d'un 
humidificateur d'air (A : vapeur d'eau, B :air). La vitesse est supposee connue, elle est orientee 
dans la direction x et ne depend que de y. La vitesse debitante (debit/section) est egale a vo. Une 
membrane separe le domaine etudie d'une ambiance ou la concentration du constituant A est 
connue et uniforme: C0 A1 (par exemple de l'air sature en humidite). La hauteur du dispositif est L. 



/\ y 



°A0 \ 


1 


4 


2 


membrane ' 









limites du domaine 



> 



air sature en humidite 
eau 



LI 



< 

a) Interieur du domaine 
conservation de la masse totale: 

dp 

+ V.(pv) = p = cste (eq. d'etat) 

at 



=> V.v = o 



3v. 



3v, 



■ + 



dx dy 



= (eq.l) v = 



v x {y[ 
o . 



repond bien a cette equation 
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conservation de la masse du constituant A: 

■ + V.(pco A v + j A ) = r A © p = cste © j A = -pD AB V(0 A (loi de transfert) © r A = 

' Ja) ^ + V.((O A v-D AB Va) A ) = 



3p(o A 
3t 



dt 



L'equation de conservation de la masse du constituant A par convection-diffusion (sans terme 
source et pour des proprietes physiques constantes) s'ecrit done : 

' :CO ^ + V.(a) A v)-D AB V 2 a) A =0 



3t 

Le terme V 2 (0 A = V.Vco A encore note Aco A est appele Laplacien de (o A • H s'agit d'une 
generalisation de la derivee seconde, caracteristique de la concavite. On obtient son expression en 
coordonnees cartesiennes en developpant le produit scalaire 5 : 



V.Vco A = 



3co A /3x 



a a 



3x ' dy) v3(0 A I dy. 



a 2 oo A a 2 co A 



3x 2 



L'expression developpee de l'equation de conservation du constituant s'ecrit done 



ao^ 
at 



■ + 



3(Q A v x aw A v 



A v y 



V 



a x 



dy ) 



D 



AB 



a co A a z co z 
"ax 1 



a y 2 ) 



= (eq. 2) 



b) Limites du domaine 

Sur la section d'entree 1 nous connaissons la fraction massique (Oao 



sur SI (x = 0, < y < L) : co, 



co 



A.O 



(eq.3) 



Sur la section de sortie 2: nous negligeons la diffusion devant la convection: 

3co, 



surS2 (x = LI, < y < L) 



3x 



= (eq.4) 



La paroi 3 est impermeable et perpendiculaire a oy: 



sur S3 (0< x< LI, y = L): 



3y 



= (eq. 5) 



Au niveau de la membrane 4, le flux convectif est negligeable et le flux diffusif est proportionnel 
a l'ecart entre les fractions massiques de part et d'autre de la membrane: 



sur S4 (0< x < LI, y = 0): p.k.(co A1 -(0 A ) = — pD AB 



3(0^ 
3y 



(eq.6) 



Si sous la membrane il y a de l'air sature en humidite a la temperature T, alors: 
x A.l =Psa t { T } / P,ot = X sa t { T } (loi d' equilibre) 



x M 

A.l A 



Ct : roA1 " x .M +(l-x ).M 

A 1 A ^ A 1 ' T 



A.l A 



A.r b 



relation entre fraction massique : co 
et fraction molaire : x 



Voir egalement le recapitulatif sur les operateurs differentiels en fin de polycopie 
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c) Conditions initiales 

II reste a definir les conditions initiales, considerons par exemple qu'initialement la fraction 
massique de vapeur d'eau dans l'air est uniforme dans le domaine et egale a celle de l'air entrant: co AO 

pour t=0 (quels que soient x et y): C0 A = C0 A0 . (eq.7) 

D'une facon generale, un probleme de mecanique des milieux continus (et plus particulierement 
le probleme de transfert de matiere pose ici ) est defini par la combinaison 

- des equations de conservations (masse totale, masse du constituant i) toujours valables 



des equations d'etat (masse volumique constante) 
des lois d'equilibre (equilibre liquide-vapeur) 



proprietes d'equilibre 
des fluides et autres materiaux 



des cinetiques de transformations (reaction chimique absente ici) proprietes hors d'equilibre 
des lois de transfert (loi de Fick) des fluides et autres materiaux 

des conditions aux limites (entree, sortie, paroi, membrane) relatives a la configuration etudiee 
et des conditions initiales (sauf en regime permanent) et aux conditions de fonctionnement 



1.6.2) Adimensionnalisation 

Definissons maintenant des variables sans dimension: 

- la concentration adimensionnelle: CO * = — - — 03 A 



03 Al -G>A.O 



elle vaut pour l'air entrant ou initialement present dans le domaine 
elle vaut 1 si l'air se sature entierement 



les coordonnees adimensionnelles: x^ 



les composantes adimensionnelles de la vitesse: v x * = — = 1, v * = — = 



v v 



En effectuant les changements de variables relatifs a cette adimensionnalisation on obtient: 
a) Interieur du domaine 



dv. * 5v v *^ 



+ ■ 



V dx* dy 



= (conservation de la masse, eq. 1 multipliee par L/vq) 



aco A * + v .L 



D AB t 
^ L 2 



D 



aco A *v x * aco A *v y *^ 



AB 



dx * 



+ 



dy : 



^3 2 co A * + 3 2 co A *^ 



dx* 2 



dy 



* 2 



= 



(conservation du constituant i, eq.2 multipliee par L 2 /D AB ) 
On voit apparaitre un temps adimensionnel appele nombre de Fourier (matiere) et note Fo. 
D AB .t t 



t* = Fo = 



L 2 L / D 



AB 



<— temps caracteristique de diffusion 



18 



Considerons un milieu semi-infini de vitesse nulle (les transferts se font uniquement par 
diffusion), la fraction massique initiale est uniforme et egale Cfl A0 . Si a l'instant t=0, on porte la 
fraction massique a la surface (x=0) a une valeur C0 A1 , alors a la distance L de la surface, la fraction 
massique reduite:((Q A -G)Ao)/(G)Ai-a>Ao) atteint une valeur de 0,48=1/2 au bout d'une duree egale a 

L 2 /D A g (demonstartion non presentee ici) 

L 2 /D AB est done un temps caracteristique de diffusion (duree necessaire pour que "l'information" 
parcoure une distance L ). 



t=0 



t=L 2 /D A 



CO 



Al 



CO 



A0_ 




infini 



On voit egalement apparaitre le rapport (v .L/D AB ) appele nombre de Peclet et note Pe. 
On peut l'interpreter comme un rapport d'ordres de grandeur entre flux convectif et diffusif. 



Pe = 



Vq- l _ P- v o(to A .i -<Q A .o) 



D 



AB 



P D 



(®A.1 -®A.o) 



flux convectif 
flux diffusif 



AB 



p. V (C0 A j — C0 A o ) represente la difference entre le flux convectif a l'entree (v ,(0 AO ) 

et celui qu'il y aurait a la sortie pour un dispositif tres long C0 A sortie — li->°° — >(0 / 

(C0 A1 -C0 A0 ) 

pD AB represente le flux diffusif a travers une epaisseur de fluide immobile L 

pour une difference de fraction massique (co A1 -co A0 ) 

L'equation de conservation de la masse du constituant i peut done 



3co A * 
3Fo 



+ Pe 



3co A *v x * + 3co A *v y 



V 



dx * 



dy : 



se reeenre: 



dx * 2 dy* 2 



J V 



= (eq. 2') 



b) Limites du domaine 

surSl (x* = 0, 0<y*<l) : co A * = 



surS2 (x* = Ll/L, 0<y*<l) 



3co A * 

dx* 
3(0. * 



= 



surS3 (0<x*<Ll/L, y* = l): —^ = 

dy * 

3co * k T 

sur S4 (0 < x* < Ll/L, y* = 0) : — = ^.(co A * -1) 

^y* d ab 



(eq.3') 
(eq. 4') 

(eq. 5') 
(rq:co A1 * = l) 
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Dans la condition aux limites relative aux transfert a travers la membrane, on voit apparaitre le 
rapport (1c.L)/Dab appele nombre de Sherwood (ou nombre de Nusselt matiere) et note Sh. Ce 
nombre compare le coefficient de transfert a travers la membrane k et D A b/L qui caracterise les 
transferts par diffusion. 

k.L k transfert a travers la membrane 

Sh = = <- 

D AB D AB /L transfert diffusif 

p.k((0 A , — C0 A ) est la densite de flux traversant la membrane pour une difference 

de fraction massique (co A1 -C0 A1 ) 

p(D AB / L). ((Q a j — (0 A o ) est la densite de flux diffusif a travers une epaisseur L de fluide 

immobile pour une difference de fraction massique (co A1 -CO A0 ) 

La condition aux limites S4 peut done encore s'ecrire: 

d(d * 

surS4 (0< x*< Ll/L, y* = 0): — = Sh((0 A *-l) (eq. 6') 

dy * 

c) Conditions initiales 

Les conditions initiales s'ecrivent simplement: 

pour t* = Fo = 0, quels que soient x * et y*: u) A * = (eq. 7') 



Pour un champ de vitesse donne, les equations 2' a 7' decrivent entierement le systeme. 

Elles permettent de calculer revolution du champ de la variable reduite Q) A * = f {x*, y*, Fo = t *} 
en fonction des parametres adimensionnels: Pe, Sh et de la geometrie: (Ll/L) 

Linteret de la forme adimensionnelle des equations est multiple: 

pour une geometrie donnee si on connait C0 A * = f{x*,y*,Fo,Sh,Pe}, on connait la solution 
pour toutes les dimensions L (Ll/L fixe), pour toutes les vitesses v , toutes les diffusivites: 
D AB , tous les coefficients de transfert: k, toutes les fractions massiques C0 A et (0 A A . 

- les resultats bibliographiques concernant des geometries simples sont faciles a manipuler. 

- on peut facilement comparer l'importance de plusieurs phenomenes; par exemple si le nombre 
de Peclet est eleve, les transferts par diffusion sont a priori negligeables devant les transferts par 
convection. 

si on cherche a caracteriser experimentalement un dispositif, le plan d'experience est beaucoup 
plus leger: il suffit de faire varier deux nombres sans dimension : Sh et Pe (en jouant par 
exemple sur le debit et la nature de la membrane) plutot que tous les parametres dimensionels: 
L, v , p, D AB , k, C0 A et (0 A1 . 

La resolution du systeme ne peut etre effectuee de fa§on analytique que dans des geometries 
simples, il faut sinon avoir recours a des methodes numeriques, de nombreux logiciels existent dans 
ce domaine. Ceci ne dispense pas de poser le probleme, par ailleurs on a toujours interet a valider la 
solution numerique dans un cas simplifie ou il existe une solution analytique. 
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Remarque: 



II est possible d'obtenir directement les parametres adimensionnels (C0A*> x *>y :| SFo,Sh,Pe) dont 
depend le probleme par une analyse de tous les parametres qui entrent en jeu 
((fl A ,(Q Asat ,CO AO ,x,y,t,L,v ,p,D AB ,k) et des dimensions fondamentales que font intervenir ces 
parametres. 

II faut remarquer, avant tout, que le probleme est lineaire (ou plutot affine) par rapport a la fraction 
massique: si C0 A * est la solution pour les deux parametres GOao=0 et C0 A i=l (toutes choses egales par 
ailleurs) , alors la solution pour des parametres CO A o et C0 A i quelconques est : 

oo A = oo AO + to A *(co A1 - ©ao). 

Ceci permet a priori de considerer uniquement la variable reduite: 00^* plutot que 

«A'«A.sat et(0 A.O 



On forme alors la matrice des dimensions des variables. 





CQA* 


X 


y 


t 


L 


v 


p 


dab 


k 


m 





1 


l 





1 


1 


-3 


2 


1 


s 











1 





-1 





-1 


-1 


Kg 




















1 





' 



Si n est le nombre de variables dimensionnelles en relation et p le rang de la matrice (dimension de 
la plus grande matrice carree extraite a determinant non nul) alors on peut exprimer la relation en 
faisant intervenir (n-p) nombres sans dimension. 

Dans notre cas n=9 et le rang de la matrice vaut 3 (colonnes L, v , p par exemple). D'ailleurs seul p 
fait intervenir une masse, il n'intervient done pas, car il ne peut etre mis en relation avec aucun 
autre parametre. 

La forme adimensionnelle de la relation fait done intervenir 6 parametres: co v A i *=f{x*,y*,Fo,Sh,Pe} 
1.7) Equation de conservation sur un domaine macroscopique (masse totale et constituants) 
1.7.1) Masse totale 



v : vitesse du fluide 



Nous allons maintenant integrer cette relation sur un systeme macroscopique dont le volume est 
note: V et la frontiere: S. Cette frontiere peut eventuellement se deplacer a une vitesse notee v s . 
On obtient evidemment : 

f — dV + f V.(pv)dV =0 y/ ^% vitesse de la frontiere 

Jv d\ Jv J 



Nous avons vu que sur un element de volume fixe on peut eenre: 
dp 

at 

accumulation bilan des flux 



+ V. (pv) 
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Deux theoremes permettent de transformer cette expression pour la rendre plus facilement 
exploitable : 

d r r 30 r 

Relation de Leibnitz : — 0. dV = — dV + v „ . n dS par exemple : = p 

(si la frontiere est fixe on peut passer la derivee sous le signe somme) 
Relation d'Ostrogradskii J (V.(p)dV= J (cp.n) dS par exemple : cp = v 



En appliquant ces theoremes on trouve : 



d_ 
dt' 



J v pdV + J s p(v-v s ). ndS =0 

accumulation bilan (sortie - entree) des flux 



Le terme J^pdV represente la masse totale contenue dans le systeme : m 

Le flux qui sort a travers une facette fixe: dS, vaut: pv.h dS . Ceci traduit bien que, si le vecteur 
densite de flux: pv , est parallele a la facette, le flux qui traverse la facette est nul et s'il est 
perpendiculaire a la facette le flux vaut: p|v|dS . ( n est le vecteur normal unitaire oriente vers l'exterieur) 




Si la frontiere est fixe, le flux de masse totale qui sort du systeme s'ecrit done : j^pv.ii) dS 
Si la frontiere est mobile, il faut prendre la vitesse relative du fluide par rapport a cette frontiere : 
v - v s . Le terme j" s p( v - v s )ii dS represente done le bilan (sortie-entree) des flux de masse dans 
le cas general (frontiere fixe ou mobile) 

Pour un systeme ayant une section d'entree et une section de sortie on peut ecrire: 



dm , . 
— + A(m) = 

dt s-e 




m : masse totale (de fluide) dans le systeme 
m e et m s : debit massique a l'entree et a la sortie 

m = l v pdV ' s A e (m) =m s -m e , m e = j^ptv - v s ).n dS 

S'il existe des apports exterieurs de masse entrant dans le systeme autrement que par les sections 
d'entree/sortie (par des parois poreuses par ex.), il faudra les ajouter comme terme source : rh ext , 

on ecrit alors : dm/dt+A(m) = rii. 
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1.7.2) Masse du constituant i 

En integrant l'equation de conservation pour le constituant i, on obtient: 
lv^T dV + LV-CP^v + l)dV = f r.dV 



J v at 

accumulation 



bilan des flux convectifs 
et diffusifs 



terme source 



En effectuant les memes transformations que pour la masse totale, on peut ecrire: 

^f v pa»,dV + ^(v-vJiidS + J s l5dS = j^dV 

accumulation bilan des flux convectifs et diffusifs terme source 
(systeme pouvant avoir des frontieres mobiles se deplacant a la vitesse v s ) 

Pour un systeme ayant une section d'entree et une section de sortie et en negligeant les flux 
diffusifs sur ces sections, on peut ecrire: 



dm 



dt 



L + A(co i .rh) = J rjdV + m^ 



mj est la masse du constituant i contenu dans le systeme 



C0 i e est la fraction massique de melange a 1' entree: co ie = 



\ se pto i ((v-v s ).n)dS 



m. 



(si la concentration est uniforme a l'entree, on a simplement: (0 ie = C0 ie ) 

m iext est l'apport exterieur de masse du constituant i (Kg/s) entrant dans le systeme 
autrement que par les sections d'entree sortie (par des parois poreuses par ex.) 

| v r ; dV represente la masse du constituant i apparaissant dans le systeme par unite de 
temps du fait de reactions chimiques 
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2) Transferts de quantite de mouvement 



Je ne sais pas ce que le monde pensera de moi. Pour ma part, j'ai I 'impression de n'avoir ete qu'un enfant qui joue 
sur la plage et se divertit en trouvant ga et la un coquillage plus joli qu'a V ordinaire, alors que le grand ocean de la 
verite reste inexplore devant moi. Isaac Newton 

La philosophic (c'est a dire la science) est ecrite dans ce tres grand livre qui se tient constamment ouvert sous nos 
yeux, I'univers, et qui ne peut se comprendre que si Von a prealablement appris a en comprendre la langue et a en 
connaitre les caracteres employes pour Vecrire. Ce livre est ecrit dans la langue mathematique, ses caracteres sont 
des triangles, des cercles et d'autres figures geometriques, sans I'intermediaire desquels il est impossible d'en 
comprendre humainement un seul mot, et sans lesquels on ne fait qu'errer vainement dans un labyrinthe obscur. 
Galileo Galilei 

2.1) Equation de conservation sur un element de volume fixe 



Dans un repere galileen, pour un point materiel, le principe fondamental de la mecanique s'ecrit 

d mv. 



pour la resultante: 



d mv 

= F 

dt 



d t 
dmv 

= F 

d t y ext 
d mv 7 
= F 

d t 



— > 

pour le moment d OMa m v 
par rapport a O d t 



— > 
= OMa F, 



d m(y.v z - 


z.v y ) 


dt 




d m(z.v x - 


x.v z ) 


dt 




d m(x.v y - 


y- v x) 


dt 



(y- F z.ext - Z - F y.ext) 
( Z -F x .ext _ X - F z.ext ) 
( X - F y.ex, -y- F x.ext) 



Ceci traduit le fait que l'accumulation (derivee par rapport au temps) de la quantite de mouvement 
selon une direction donnee est egale a la composante des forces exterieures selon cette direction. 
Ces forces exterieures apparaissent done comme une source de quantite de mouvement. 

Au paragraphe 1.1.2, nous avions ecrit le bilan de la masse du constituant i sur un element de 
volume fixe: 

3p(0 ; 



3t 

accumulation 



+ V^pcOjV) + 

bilan des flux 
convectifs 



V.Q) 

bilan des flux 
diffusifs 



terme source 



La fraction massique C0{ , qui etait la variable transported, representait la masse du constituant i 
par unite de masse totale. 



La quantite transportee a laquelle nous nous interessons maintenant est la quantite de mouvement 
dans une direction donnee (ox: par exemple). La variable transportee correspondante est 
simplement v x qui represente la quantite de mouvement dans la direction x par unite de masse 
(totale). 
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En remplacant (0[ par v x , on peut ecrire les differents termes du bilan de quantite de mouvement 
dans la direction x sur un element de volume fixe: 



3t 

accumulation 



+ 



bilan des flux 



V. (pv x v) = 



forces volumiques 
a distance 



+ 



par les elements de fluide 
voisins a travers les facettes 



bilan des forces exercees 



convectifs 



est l'accumulation de quantite de mouvement (dans la direction x) par unite de 
temps (divise par le volume de l'element) 



V.(pv x v) represente le bilan des flux convectifs de quantite de mouvement (dans la direction 
x) qui sortent de l'element par ses 6 facettes. 



Le fluide qui traverse une des facettes avec une vitesse v transporte avec lui (par 
unite de volume) une quantite de mouvement (dans la direction x) egale a p.v x ; de 
la meme facon qu'il transporte une masse de constituant i egale a p(Qj 

Le terme pv x v est done le vecteur densite de flux de quantite de mouvement dans 
la direction x; de la meme facon que pWjV est le vecteur densite de flux de masse 
du constituant i. 



La force volumique la plus couramment rencontree est la gravite. Dans ce cas 
g = [g x ,g y ,g z ) avec g x =g y =0etg z = -g : acceleration de la pesanteur. 



II reste a ecrire le bilan des forces exercees (dans la direction x) par les elements de fluide voisins 
a travers les facettes. Raisonnons en deux dimensions. 

La force qui s'exerce a travers une facette depend de l'orientation de cette facette (et de la position 
de cette facette). On appelle contrainte totale, la force par unite de surface (homogene a une 
pression). On adopte les conventions suivantes : 

71 xx composante selon x de la contrainte exercee a travers une facette perpendiculaire a ox. 
7C xy composante selon y de la contrainte exercee a travers une facette perpendiculaire a ox. 

(dans le sens des x croissants, contrainte exercee par A sur B sur la figure ci-dessous) 



PSx 



represente les forces exterieures a distance par unite de volume. 



y 




X 



A 



B 
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7t yx composante selon x de la contrainte exercee a travers une facette perpendiculaire a oy 
7t yy composante selon y de la contrainte exercee a travers une facette perpendiculaire a oy. 

(dans le sens des y croissants, contrainte exercee par A sur B sur la figure ci-dessous) 









7Cyy 






B 







yx ' yy- 



7ly x 



Faisons maintenant le bilan des forces exercees dans la direction x sur le fluide contenu dans un 
element de volume a travers ses facettes. 

-7t yx {x, y+Ay/2} 



7l xx {x-Ax/2 , y} 



+ 



- 7C xx {x+Ax/2 , y} 



7C yx {x, y-Ay/2} 

La force exercee par le fluide situe a gauche sur celui de l'element vaut : Ay.Az .7i xx {x-Ax/2 , y} 
(composante selon x, facette _L a ox, sens des x croissants) 

La force exercee par le fluide situe a droite sur celui de l'element vaut : - Ay.Az .7l xx {x+Ax/2 , y} 
(sens des x decroissants, on applique la loi de Taction et de la reaction , d'ou le signe negatif) 

La force exercee par le fluide situe en dessous sur celui de l'element vaut : Ax.Az .TC yx {x , y-Ay/2} 
(composante selon x, facette _L a oy, sens des y croissants) 

La force exercee par le fluide situe au dessus sur celui de l'element vaut : - Ax.Az .7U yx {x , y+Ay/2} 
(sens des y decroissants — > signe negatif) 

Faisons la somme et divisons par le volume de l'element (Ax. Ay.Az) on obtient : 

7l xx {x-Ax/2,y}-7l xx {x + Ax/2,y,z,t} 7l yx {x,y - Ay / 2} - 7l yx {x,y + Ay} 



Ax. 



Soit en passant a la limite : - 



■ + ■ 



Ay 



■ + ■ 



3x 8y 
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Sous forme developpee en deux dimensions 1' equation de conservation de la quantite de 
mouvement dans la direction x s'ecrit done : 



3pv x _ ( 3pv x v x _ dpw x w y ) ( d %xx d% yx 



at 



■+ 



En notant 7t v = 



3x 



+ 



+ 



+ 



dy J v 3x dy J 



Pg> 



on peut contracter 



+ - yx sous la forme V.(S v ) 



3x 3y 



D'ou la forme contracte de l'equation de conservation de la resultante de quantite de mouvement 
selon ox : 



3t 

accumulation 



+ V. (pv x v) + V. (fij 

bilan des flux bilan des forces 
convectifs surfaciques 



Pgx 

forces volumiques 
a distance 



V.(tc >x ) represente le bilan des forces qu'exerce l'element de volume sur ses voisins au 
travers de ses 6 facettes dans la direction x (divise par le volume de l'element). 

Cela correspond a un bilan (sortie - entree) de flux de quantite de mouvement 

Des equations equivalentes vont bien sur exprimer la conservation de la quantite de mouvement 
dans les direction oy et oz 

Une analyse de l'equilibre en rotation d'un petit element de volume autour de l'axe oz aboutit a 
7ly X =7t xy (non demontre ici). L'equation de conservation du moment de la quantite de mouvement 
se traduit done par le fait que l'ordre des indices est indiferent pour la contrainte totale: 



Ttyx — 7C X y , TCyz — T^zy , ^zx — ^xz 



Les differentes composantes des contraintes peuvent etre regroupees dans un matrice symetrique 
appellee tenseur des contraintes totales : 





Xx 






71 = 


^yx 


% yy 






V^zx 







Attention: 

En coordonnees cylindriques, les vecteurs de base dependent de la position du point, ce qui rend 
l'ecriture des operateurs gradient et divergence plus complexe. II ne suffit done pas de remplacer 
x par r et y par pour passer en cylindrique. On trouvera 1' expression des equations de bilan dans 
differents systemes de coordonnees par exemple dans Transport phenomena (cite en preambule). 

De tres nombreux auteurs utilisent une convention de signe differente pour les contraintes 
totales, ils notent ay = -Try. La convention adoptee dans ce cours est celle utilisee dans Transport 
Phenomena, elle a l'avantage de pouvoir assimiler les forces a des flux de quantite de mouvement 
et de conduire a des lois de transfert similaires pour la quantite de mouvement (fluide newtonien 
incompressible), la matiere (loi de Fick) et la chaleur (loi de Fourier). 
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Ecrivons maintenant la conservation de la quantite de mouvement dans les trois directions x,y et z. 



dpv x , ( 3pv x v x apv x v y a P v x v z ^ f 3tc xx 3ti yx 3ji 



3x 



3y 



+ 



3z 



3p\ 



at 



■+ 



apv y v x apv y v y apv y v 



dy 



+ ■ 



dz 



+ 



+ 



■ + 



+ ■ 



dx dy dz 



■ + ■ 



3tc„ 



dx dy 



+ 



3pv z _ ( 3pv z v x _ apv z v y a P v z v z ^ _ ( dn xz _ 3ti 



at 



+ 



3x 



3y 



3z 



+ 



dz 
5tc 



+ ■ 



dx dy 



■ + ■ 



3z 



= PSx 
= Pgy 
Pgz 



Ceci peut s'ecrire matriciellement: 



^(pv x , pv y , pv z ) + 



3x ' 3y ' dz 



f Pv x v x pv x v y pv x v\ 

P V y V X P V y V y P V y V z 
\P V z V x P V z V * P V z V z/ 



+ 

dx dy dz 

On voit apparaitre deux matrices: 

r \ 
pv v pv v pv v 

X X X y X ; 







TC 






TC 






d_) 


v^ zx 


TC 



xy 



TC 

yy yz 



zy 



gx> Pgy. Pgz 



pv v pv v pv V 

y x y y y z 

pv v pv v pv V 

Vz X zy z z J 



notee pvv 



f 






TC 


TC 


TC 


XX 


xy 




TC 


TC 


TC 


y X 


yy 




TC 


TC 


TC 


V zx 


zy 





notee TC 



zz / 



Rappelons la forme du vecteur symbolique nabla: V = 



3 3 3 
V3x 3y 3z 



(coord, cartesiennes) 



Les trois equations de conservation de la resultante de la quantite de mouvement peuvent done etre 
resumees sous la forme 



dpv 
accumulation 



+ V. (pvv) 

bilan des flux 
convectifs 



+ V.(TC) 

bilan des forces 
surfaciques 



Pg 

forces volumiques 
a distance 



Les trois equations de conservation du moment de la quantite de mouvement s'expriment, quant a 
elles par le fait que la matrice TC est symetrique: 



*xy 


= TC 




^yz 


= TC 




^z X 


= TC 



<=> TC = (tc)' 



ou (Tc)'est la matrice transposed de TC obtenue en 
inversant les indices des lignes et des colonnes 



La matrice TC est appelee tenseur des contraintes totales. Dans un repere cartesien on peut 
l'assimiler a une matrice symetrique 
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2.2) Equation de conservation sur un element de masse fixe suivi dans son mouvement . 



Effectuons la meme operation qu'au paragraphe 1.2, c'est a dire soustrayons l'equation de 
conservation de la masse multiplied par v x a celle de conservation de la quantite de mouvement 
dans la direction x: 

' 3p + 3pv i+ aPX L+ 3pv t = @ (_ v ) 



■ + ■ 



dt 3x 3y dz 

© 



apv x apv x v x apv x v y a pVx v 



at 



-+ 



•+ 



X z 



3y 



a z 



37t_ a^ 3tc, 



d'ou 



av„ av v av v av 

+ v. — - + v„ — - + V 



at 



a x 



soit vectoriellement : p : 



Dv. 



y a y z a z 

dv. 



ax ay az 
+pg 



+pg x 



Dt 



at 



-+(v.vK 



a7t xx | a7t yx | d% z 

dx dy dz 
-V.7t x +pg x 



en coordonnees cartesiennes et sur 1' ensemble des trois composantes de la vitesse : 

"av 



Dv 

P = P 

Dt 



at 



■+ 



(v.v)v 



= -V.7C + pg 



La derivee particulate de la vitesse Dv / Dt (evolution de v en suivant le mouvement) n'est autre 
que 1' acceleration d'une particule materielle de fluide. Le terme de droite de l'equation represente 
les forces exercees sur cette particule a travers ses facettes et les forces a distance (par unite de 
volume). 



bilan sur un element de volume 
fixe dans l'espace 

systeme ouvert : du fluide entre et sort a 
travers les facettes et transporte avec lui de 
la quantite de mouvement (par convection) 



apv 

"aT 



+ V.(pvv) = -V.7l + pg 



bilan des flux 



accumulation + (convectifs) a = 
travers les facettes 



terme source : 
forces exterieures 



bilan sur une particule de fluide 
suivie dans son mouvement 

systeme ferme : la particule se deplace 
et se deforme avec le fluide, il n'y a pas 
d'echange de flux convectif par les facettes 



Dv 
Dt 



= - V.71 + pg 



masse * acceleration = forces exterieures 
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Remarque : la derivee particulaire est un cas particulier de derivee d'une fonction de plusieurs 
variables en suivant un chemin . 

exemple : y=f{x,t} 

si Ton suit un chemin caracterise par x=g{t} alors y=f{g{t},t} on peut ainsi definir la derivee 
de y par rapport a t en suivant ce chemin : 



dy 
dt 



d[f{g{t},t}] 



=g{t} 



dt 



on peut egalement raisonner en differentiant y selon un chemin ou dx=g' {t}dt 



dy dy i 
dy = T L dx + ^ L dt -» dy I ,„ 
: dx dt • 7 idx= g '.dt 



— > 



dy 
dt 



dx=g'.dt 



ax~ g{t} + ¥ 



cas de la derivee particulaire de v: v = v{x,y,t} chemin dx=v x .dt et dy=v y .dt 

dv 



Dv 
Dt 



dv 
dt 



dx=v x .dt 
dy=v y .dt 



dx x 



dv dv 
dy~ Vy+ ^t 



En coordonnees non cartesiennes 1' expression de la derivee particulaire ci - dessus n'est plus 
correcte du fait que les vecteurs de base ne sont pas fixes, 1' expression generalisee s'ecrit: 

Dv _ 

at 



P 5T = P 



+ — V(v.v) + (v a v) a v 

^ 2 v ; 



= -V.Ti + pg 



oil v a w represente le produit vectoriel de v et de w 

et V a v represente le rotationnel de v qui en coord, cartesiennes vaut: 



3v z 




3y 


dz 






dz 


dx 






dx 


3y, 



2.3) Equation d'etat 

Nous avons deja indique que la masse volumique est, en general, fonction de la pression et de la 
temperature, citons quelques cas particuliers. En premiere approximation, la masse volumique des 
liquides peut etre consideree comme constante. II en est de meme pour les gaz, si la pression et la 
temperature varient peu et que les vitesses sont faibles devant la celerite du son, c'est le cas pour 
beaucoup de problemes d'ecoulement d'air (v<60m/s). 

Pour un gaz parfait de masse molaire M, nous connaissons une relation qui lie p,p et T: 
pV = nRT => p = RpT/M (p = nM/V) 

Si la temperature est constante (gaz parfait isotherme) on obtient une relation immediate entre la 
pression et la temperature: : p / p = RT / M = cste 



Dans certains cas on peut considerer que le gaz parfait s'ecoule sans echanger de chaleur et sans 
frottement (gaz parfait isentropique), on sait qu'alors pV Y = cste, d'ou: p / p 7 = cste 
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2.4) Loi de transfert 6 

2.4.1) Fluide au repos - Contrainte de pression 

Un fluide est defini par le fait qu'au repos, les contraintes sur une surface elementaire sont 
normales (perpendiculaires) a cette surface et sont independantes de l'orientation de la surface. La 
contrainte (force par unite de surface) est alors egale a la pression du fluide en ce point. 



contrainte exercee 
B par A sur B 



B -k- 



71 yy =p 



71 XX =p 



A contrainte exercee par A sur B 



Ainsi pour un fluide au repos : 7t xx = 7t yy = 7C ZZ = p et 7C xy = % y 



71, 









*xx 


11 xy 


^xz 




( 

P 





& 




f \ 





°1 


Soit: 




71 yy 


% yz 







P 





= P 





1 







V^zx 


^zy 






vO 





PV 




vO 





h 



en condense : tc = p I 



Si on fait le Man des forces dans la direction x sur un volume elementaire, on voit que les forces 
de pression rapportees a l'unite de volume valent: 

(p{x-Ax/2, y, z}-p{x + Ax/2, y, z}).Ay.Az 3p 

— - — soit en passant a la limite: 

Ax.Ay.Az 3x 



\1/ 



p{x- Ax/2, y, z} 




p{x+A x/2, y, z} 



En generalisant aux trois composantes, on voit que la resultante des forces de pression, rapportee a 
l'unite de volume, est l'oppose du gradient de la pression: 



dp/dy 
dp/dzj 



-Vp 



C'est bien ce que Ton obtient en prenant l'oppose de la divergence de 7C (fluide au repos) 

-V.I = -V.pl = -Vp 



En mecanique, on parle plutot de loi de comportement du materiau 
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Equilibre statique d'un fluide 

Pour un fluide au repos et soumis a la gravite, le bilan de quantite de mouvement s'ecrit: 











V.7i = pg <=> V.pl = pg o 


dp/dy 


= P 







v3p/3zy 




v-gJ 



Si p est constant : p = p{z = 0} - pgz = p{h = 0} + pgh (h = -z) 
2.4.2) Fluide en ecoulement - Contraintes de frottement visqueux 

Lorsque le fluide est en mouvement, en plus des contraintes de pression, il s'exerce des contraintes 
de frottement visqueux notees x xx , x yx , x zx ,x xy etc. Ces contraintes sont regroupees dans une 

matrice (symetrique) appelee tenseur des contraintes de cisaillement (on parle egalement de 
contraintes visqueuses ou de contraintes de deformation). 



f 


^xy 


^xz 




c 

p 





0^ 




XX 


X 

xy 


^xz 




71 yx 


TC yy 


^yz 







p 





+ 


X 

yx 


X 

yy 




en condense 71 = p. I + x 




^zx 






v0 





Pv 








X zzV 





contraintes contraintes contraintes de cisaillement 

totales de pression (ou contraintes visqueuses) 



Soit, par exemple : 7t xx = p + x xx ou 7t yx = x xy 
Ecoulement unidirectionnel - Cisaillement simple 

Considerons le cas du cisaillement simple. Du fluide est contenu dans l'entrefer de deux plaques. 
Celle du bas est fixe, celle du haut est en translation selon x a une vitesse v max . 



y 















e 


X ' %x 





Alors la vitesse varie lineairement selon y 



e dy e 

dv 

— - est appele taux ou vitesse de cisaillement 

dy 
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Une couche de fluide exerce sur celle situe au-dessus (et la couche de fluide superieure exerce sur 
la plaque du haut) une force de frottement visqueux dans la direction opposee a l'ecoulement. 
Cette force par unite de surface est la contrainte de cisaillement (ou de frottement visqueux) : x yx 

(la force est exercee a travers une section perpendiculaire a Oy dans le sens des y positifs et elle 
est dirigee dans la direction Ox). 

dv 

Pour un fluide newtonien on a : x v = -pv — - 

yx dy 

v est appele viscosite cinematique (m 2 /s) 
(i=pv est appele viscosite dynamique (Pa.s) 

Noter l'analogie avec les transferts de matiere (loi de Fick) et de chaleur (loi de Fourier) 

^^^^^rarjide^^^^^^ ^ ^ ^^^^ojicentr^^^^^ 

Xy x flux de ^qte de mouvement j q y flux^de chaleur jA.y flux ^de matiere 

lent ^^^^^rad^^^^^^™ ^^^^^^Tlue^^^^^" 



dy 



J 



--I — 

l| - y dy 



dH 

pa." — (C = cste) 

dy p 



j Ay = -PD 



AB 



dco^ 
dy 



v , a = X I (p.C p ) et D AB ont la meme unite : m 2 /s 

On peut assimiler la contrainte de cisaillement a un flux diffusif de quantite de mouvement et la 
viscosite cinematique a un coefficient de diffusion de la quantite de mouvement . Du fait de 
1' agitation moleculaire (schema ci-dessous) des particules vont traverser le plan perpendiculaire a 
oy. Autant de particules iront de haut en bas que de has en haut (conservation de la masse). Mais 
les particules qui descendent transportent avec elles une quantite de mouvement (selon ox) plus 
importante que celles qui montent, d'ou un flux net de quantite de mouvement. 

Les particules rapides qui descendent vont entrer en collision avec les particules plus lentes de la 
couche inferieure et reciproquement pour les particules qui montent, ceci explique la force 
visqueuse d'entrainement de la couche lente par la couche rapide. 



transfert diffusif de quantite de mouvement 

A 



transfert diffusif de matiere 



o o 



o o 



o v o 



D'ailleurs le nombre de Prandtl: Pr= v/a est proche de 1 pour l'air (Pr{air}=0,71). 
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La viscosite dynamique est fonction de la temperature: elle augmente avec la temperature pour les 
gaz et diminue pour les liquides, mais elle varie peu avec la pression des gaz. 



Attention de nombreux fluides rencontres dans le domaine biologique, agro-alimentaire et 
environnemental sont non-newtoniens (compote, ketchup, boues, sang a l'echelle des capillaires) 

2.4.3) Contraintes dans un fluide en ecoulement quelconque 



Ce que nous venons de voir n'est valable que si l'ecoulement est unidirectionnel: v y =v z =0 et que v x 
ne depend que de y. 

Dans le cas des transferts de chaleur ou de matiere, si T (resp.G)) est une fonction quelconque de x, 
y et z. On generalise la relation en utilisant la notion de gradient (pente generalisee). 

Dans le cas d'un ecoulement quelconque, on serait tente de faire de meme pour les transferts de 
quantite de mouvement. Cela conduirait a : 



dv x dv x dv 

txx =-pV^ > V ="PVV L » X zx =-p V 



dx 



dy 



dz 



txy ="PV- 



dx 



etc. 



(Faux) 



Nous savons que la matrice x est symetrique , on doit avoir par exemple x yx = x xy 

dv 3v 
Or il n'y a aucune raison a priori pour que — — soit egal a — — . 

dy dx 

On ne peut done pas simplement raisonner par analogie, il faut revenir a la cause premiere des 
contraintes de frottement visqueux : la deformation des particules de fluide. Si le fluide se translate 
en bloc ou tourne en bloc, comme un solide, le fluide ne se deforme pas et il ne s'ajoute pas de 
contraintes supplementaires a celles de pression. On decompose done localement le mouvement en 
une combinaison de translation, de rotation et de deformation. Puis on enonce que seule la 
deformation contribue aux frottements visqueux et on considere une dependance lineaire pour les 
fluides newtoniens. 

Finalement pour un fluide newtonien incompressible , on a : 



^xy = V = "PV 



yx 



dv„ dv 



- + ■ 



dy dx 



t X x =-2pv 



dx 



et des relations similaires pour x yz =x zy , x zx = x x 



Xyy ^ X, 



De fa§on mnemotechnique, on peut retenir qu'on 'symetrise' la relation pour x xy =x yx et que pour 
x xx on prend la meme relation que pour x xy en changeant l'indice y en x, d'ou le facteur 2 

En deux dimension on peut ecrire matriciellement : 



r 



x 

V y x 



= -pv 



3v x 
2 — - 



3v x 3v 
— - + — - 

dy dx 

3v„ 



dy dx 



dy 



rappel (J, = pv 
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Le paragraphe suivant presente cette demarche en 2 dimensions, 
instant t instant t+dt 




translation 



rotation 



^ — ^ deformation 

Considerons une particule de fluide qui a l'instant t forme un carre de cote dl 
A : coin inferieur gauche est situe en ( x , y ) 
B : coin inferieur droit est situe en ( x+dl , y ) 
C : coin superieur gauche est situe en ( x , y+dl ) 

A l'instant (t+dt) cette particule de fluide s'est deplacee, a tourne et s'est deformee 
Un point M de coordonnees (x M , Ym ) quelconque de la particule de fluide 
devient M' de coordonnees ( x M +v x {x M ,yivi}dt , y M + v y {x M ,yivi}dt ) 

A devient A' qui est situe en (x +v x {x,y}dt , y + v y {x,y}dt ) 
B devient B' qui est situe en ( x +dl + v x {x+dl,y}dt , y + v y {x+dl,y}dt ) 
C devient C qui est situe en (x + v x {x,y+dl}dt , y +dl + v y {x,y+dl}dt ) 




coord, de C = 



y + dl + 



v x {x,y}+dl-^4x,y}Jdt 
v y {x,y} + dl^{x,y}jdt 



Les vecteurs qui determinent les cotes 



sont a l'instant t 



AB = dl 



AC= dl. 



et deviennent a l'instant (t+dt) : A' B' = dl. 



< dw x > 
1 + ^.dt 
dx 

-r-^.dt 
^ dx V 



dt 



A'C'=dl. 



dy • 

3v. 



1 + 



dy 



.dt 



La matrice de changement de base s'ecrit done : 

r dw x 3v x ^ (7n, dv ^ 

1 + — ^.dt -^H-.dt 



3x 



3x 



3v, 



.dt 1 + 



dy 



At 



= 1 + 



dx dy 



3x 3y 



dt 



En compos ant Vet v on peut former une matrice 

fdv. dv > 

'd/dx ^ 



Vv 



d/dy 



(v x ,v y ) 



3x 3x 
dv„ 3v v 



la premiere colonne est Vv x 
la seconde colonne est Vv, 



dy dy 

La transpose de Vv est appelee gradient du vecteur vitesse (Vv) 1 = 
La matrice de changement de base s'ecrit done : I + (Vv) .dt 

La matrice de changement de base d'une rotation suivi d'une deformation lineaire symetrique 
s'ecrit: 

A ^cosCa)) - sin((0) A 



5v x 


3v x ^ 




a y 


av y 




dx 


ay j 



'l + a c 
V C 1 + by 



Vsin((0) cos((0) 



Si la deformation et la rotation sont de faible amplitude (a, b, c et (0 petits), on a 
cos((0) ~ (l-(0 2 12) ~ 1 ,sin(co) ~ 0), d'ou : 



1 + a c Y 00 ^* ) -sin((0)^j 
c l + bAsin((Q) cos((0)J 



1 + 



a c 
c b 



1 + 



-(0 
(0 





U c^ 




'0 




I + 




+ 




J 




VC by 




V(0 



((l + 8l)(l + 82) = 1 + 81 + 82) 



symetrique antisymetrique 



La matrice de changement de base decrivant la rotation et la deformation de la particule de fluide 
s'ecrit : I + (Vv)\dt . 
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D'ou l'idee de decomposer (Vv) en une partie symetrique plus une partie antisymetrique 



3v y 3v y 



dx 3y 



dt = 



1 



( 3v s 



3v 3v x 

— - + — - 
V dx dy 



3y 3x 
oy y 



l 





dV y 9V X 

I, 3x 3y 



a Vy ^ 



3y 3x 
.0 



f a c\ (0 -0^ 



\c bj 



\® J 



(Vv) 1 



1 



-(Vv + CVv) 1 ) 



1 - f - 
+ t((Vv)'-Vv) 



= d + a 



On peut done identifier: 

- le vecteur de translation 



la matrice liee a la rotation 



— (0 
0) 



a: v{x}dt 

r 

1 



a: 







dv dv v 



3v y av x 



la matrice liee a la deformation 



1 



a c 
c b 



a : 



dx dy 
dx 

dv„ dv v 



dy dx 
.0 



dt = Q dt 



dv dv 
+ 



+ - 



dy dx 



dx dy 



dt= Ddt 



La matrice D = ~(Vv + (Vv) ) , partie symetrique du gradient de vitesse, est appelee tenseur des 

taux de deformation (ou des vitesses de deformations) e'est de lui seul que dependent les 
contraintes de cisaillement qui s'ajoutent aux contraintes de pression. 

La relation la plus simple que Ton puisse imaginer entre les contraintes de cisaillement et les taux 
de deformation est une relation lineaire. On retrouve ainsi sous une forme condensee les equations 
qui caracterisent les contraintes de cisaillement (ou visqueuses) des fluides incompressibles 
newtoniens. Le coefficient de proportionnalite est : -2 |i 



x = -2|iD = -2(1 



1 



-(Vv + (Vv) 1 ) 



c.a.d. 



fx x 

XX 

V X yx 



X X / 



dv v dv, 

- + - 



dv. 



+ ■ 



dy dx 



dx dy 

3v y 

dy 



Verifions que le tenseur des taux de deformation D est nul pour une translation en bloc ou une 
rotation en bloc et que Ton retrouve le resultat precedent pour un cisaillement simple. 



pour une translation en bloc on a: 



v x = a = cste 
v = (3 = cste 



d'ou (Vv)' = 



dv dv 



dx dy 
dv y dv y 



f° °1 
v0 0y 



dx dy 

le gradient de vitesse est nul et done a fortiori D 
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pour une rotation en bloc autour de (xQ,yo) on a: 



v x =-©(y-y ) 



d'ou (Vv)' = 



3v. 



dx dy 
dv y dv y 



-w 



V© 



v y = to(x-x ) 

dx dy 

le gradient de vitesse est antisymetrique, D est done bien nul 
pour un cisaillement simple on a: 



dv 

v x =a-y = — -.y 



dy d'ou Vv 



v y =0 



f o 


0^ 




'o dv 

dy 


dv x 





(Vv) 1 = 


I dy 


J 




v J 



->Vv + (Vv)' = 





dy 



dv^ 
dy 





loi de Newton x = -|i(Vv + (Vv) 1 ) -> 



x x 

xx xy 
T T 

V y yy J 





dv. 



dy 



dv^ 

dy 
o 



On retrouve bien la contrainte de cisaillement x yx = -|x(dv x / dy) 

Les calculs sont plus lourds mais sont generali sables sans difficultes a trois dimensions, dans ce 
cas, le tenseur des taux de deformation vaut: 



D = -(Vv + (Vv)')=- 
2 2 



2^ 
3x 

dv^ 3v N 



dv v dv, dv, dv. 



■ + - 



■ + - 



9x 3y 3x dz 



3v 7 3v y 

■ + - 2 - + - 

dy dx dy dy dz 

3v„ 3v, 3v, „3v. 



+ - 



+ - 



3z 3x 3z 3y dz 

Finalement, pour un fluide newtonien incompressible, on obtient l'expression ci-dessous pour le 
tenseur des contraintes totales. 

7C = pi + x avec x = -2(iD = -2(i 



I(Vv + (Vv) 1 ) 



1 : tenseur des contraintes totales 
p : pression statique 
X : tenseur des contraintes de deformation 





^xy 






f v 





0^ 






X xy 


Xxzl 


^yx 


*yy 


*yz 







p 





+ 


X y* 


X yy 


X yz 


v ^zx 


^zy 


^zzy 




v° 









X 

^ zx 


X zy 


X 

zz J 



avec 



yx 



'xy 



yy 



zy 



yz 



3v x 

2 — - 
dx 

3v„ 3v v 



3y 
3z 3x 



dx 
3v. 



3v v dv. 



3x 3y 
3v„ 



3y 

— - + — - 
dz dy 



3v 3v 
— - + — - 
dx dz 

3v, 3v v 



+ 

dy dz 

o 3V z 

2 — - 
dz 



Si Ton voulait traiter des fluides compressibles, il faudrait encore decomposer la deformation en 
une expansion uniforme dans toutes les directions et une deformation a volume constant. 
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2.5) Conditions aux limites 



a) Conditions d' entree et de sortie 



Suivant le probleme etudie on peut ecrire des conditions aux limites portant sur la vitesse ou sur la 
pression. 

Conditions portant sur la vitesse : 
Sur la section d'entree le profil de vitesse du fluide peut etre connu. 

Sur la section de sortie la vitesse est rarement connue (dans ce cas il faut veiller a ce que la 
conservation de la masse soit verifiee globalement entre l'entree et la sortie) 

Si la sortie est placee sur une conduite de section constante et suffisamment loin du debut de 
cette conduite, on peut considerer que l'ecoulement est etabli. Cela signifie que le profil des 
vitesses ne varie plus dans le sens de l'ecoulement. Par exemple, s'il s'agit d'une conduite dont 
l'axe est oriente selon ox (section de sortie perpendiculaire a ox) on pourra se contenter d' ecrire : 



x y_ z q 

dx dx dx 

On peut montrer que dans ce cas la vitesse est orientee selon l'axe de la conduite et que la 
pression est uniforme sur la section. 



Conditions portant sur la pression 



II existe des cas ou Ton connait plutot la pression sur la section d'entree et/ou de sortie. Lorsque 
Ton a affaire a un reservoir de liquide en contact avec un gaz, par exemple, il est commode de 
placer l'entree (ou la sortie) juste en dessous de l'interface liquide/gaz ou la pression est connue. 
Si le reservoir est en contact avec 1' atmosphere, la pression a la surface libre est connue, elle est 
egale a la pression atmospherique. 



b) Conditions de paroi 



Au niveau d'une paroi fixe, pour un fluide visqueux, le vecteur vitesse est nul (condition 
d'adherence). Non seulement la composante normale a la paroi, mais egalement les composantes 
tangentielles sont nulles. Si la paroi est mobile, la vitesse du fluide est egale a celle de la paroi. 

Remarque : dans le cas de 1' approximation de fluide parfait (viscosite nulle), le fluide peut glisser 
sur la paroi, on ecrit alors seulement que la vitesse normale a la paroi est nulle : v.n = . 
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2.6) Combinaison des equations-Adimensionnalisation 



Nous allons maintenant combiner les differentes equations dans le cas d'un fluide incompressible ( 
p=cste) newtonien de viscosite constante (|i=cste) et de composition uniforme. Prenons 
l'exemple schematise ci-dessous et raisonnons en deux dimensions. 



> x 



8 



L2 



limites du domaine 



10 



L3 



L1 



2.6.1) Combinaison des equations 
a) Interieur du domaine 



conservation de la masse © loi d'etat 

^ + V.(pv) = © p = cste 



=> V.v = o 



3v T dv v 



■ + 



3x 3y 



= (eq.l) 



conservation de la quantite de mouvement © loi de transfert : 

Commencons par exprimer le tenseur des contraintes totales (pression+deformation) dont nous 
allons avoir besoin dans 1' equation de conservation de la quantite de mouvement. 



fluide visqueux : 71 = pi + x 



fluide newtonien : 



y yx yy J 



+ 



x x 

xx xy 
X X 

V yx yyy 



't X ^ 

xx xy 

X X 

V y x yy V 



-V- 



2 — - 
3v„ dv, 



3y 3x 



dv. 3v. 



3y 3x 

2 — L 
3y 



d'ou 



-2|x 



dv. 



3x 
3v, 



3x 



3v„ dv x 



- + - 



^xy =^yx 



3y 3x 
3v 



3v. A 



- + - 



3y 3x 



o 3v y 

Sy=-^" 



dy 



Kyy =P-2|I- 



^y 

3y 
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Nous pouvons exprimer 1' equation de conservation de la quantite de mouvement, soit composante 
par composante, soit vectoriellement . Commencons par la premiere methode. 

Conservation de la composante selon ox de la quantite de mouvement : 

d(pv x ) a, - _ , 3(pv x ) d(pv x ) ^(P v x v y) dn xx d7t 

+ V.( P V x V + 7U. x ) = pg x O + + _^_ + _iHL + _2_ = pg x 



at M x - x/ ,ox at 

Remplacement de jc xx et jc yx par leurs expressions 



dx 



3y 8x 3y 



d(pv x ) d(pw 2 x ) 3(pv x v y ) 3 p ( 3 2 v x d 2 v x d 2 v y ^ 



3t 



- + 



3x 



■ + 



dy 



3x 2 dy 2 3y3xy 



+ Pg x 



dp 



a 2 v. d 2 v a ( a v , 3v v v * 



■+--—+ 



dx 2 dy 2 dxv 3x 3y 



■ + ■ 



JJ 



+ Pg> 



d(pv x ) d(pv 2 x ) 5(pv x v y ) dp 



3t 



■ + 



3x 



■ + ■ 



dy 



dx 



d 2 v d 2 v. 



+ 



3x 2 dy 2 y 



conservation de la masse 
' v~ ■ 





+ Pg x (eq-2) 



En utilisant le vecteur nabla cette equation s'ecrit: * + V.(pv x v) = -— + pvV 2 v x + pg x 

dt dx 

Ecrivons maintenant directement la conservation du vecteur quantite de mouvement et combinons 
la avec la loi de transfert (cf. tableau synthetique des relations sur les operateurs en fin de 
polycopie). 

dp\ ■ = - 

-^- + V.(pvv + 7t) = pg ©p = cste© 7U = pI-pv(Vv+(Vv) t ) 

=> ^ + V. (pvv) = -V. (pi) + pvV. ( Vv + ( Vv) 1 ) + pg 

or V.(pI) = Vp , V.(Vv)' =[ V(V.v) ]' =0 et V.(Vv) = V 2 v = (v 2 v x , V 2 v y ) 
D'ou la forme condensee des equations des 2 composantes de la quantite de mouvement: 



d pv - 

— — + V. ( p vv) = - Vp + pv V 2 v + pg <=> 
at 



3pv x - dp 

+ V.(pv x v) = + pvV 2 v x + pg x 

dpv y - dp 

— + V.(pv y v) = -— + P W 2 v y + pg y 



d 2 v x d\ 
dx 2 + dy 2 

(generalisation en plusieurs dimensions de la derivee seconde caracteristique de la concavite d'une courbe) 



Le terme V 2 v v = ^- H ^- encore note Av v represente le Laplacien de v x 

x dx 2 dy 2 x 
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Notez l'analogie avec l'equation de conservation de la masse d'un constituant 



^^ + V.(pa) i v) = pD i V2(o i +r i <-> + V.(pv x v) = pvV 2 v x -^ + pg x 

at at dx 



masse d' un constituant (cf. §1.6) 



quantite de mouvement selon ox 



Finalement, pour un fluide incompressible, newtonien de viscosite constante, les equations 
de conservation de la masse et de la quantite de mouvement, appelees equations de Navier- 
Stokes s'ecrivent : 

Forme conservative developpee: 



3v x 5v y n 

— -+ — - = 

dx dy 










3(pv x ) | 3(pv 2 x ) | 3(pv x v y ) 
at dx dy 


dp 
~dx- + » 


v dx 2 


o 2 v x ^ 

+ — - 

dy 2 J 


+ Pg x 


3(pv ) 3(pv x v ) o(pv^) 


dp 
"V* 1 




d 2 Y ^ 

+ y 


+ Pgy 


+ + = 

dt dx dy 


v ax 2 


dy 2 j 



Forme conservative condensee 



V.v = 



+ V.(pvv) = -Vp + uV 2 v + pg 

dt 



Forme convective developpee (obtenue en retranchant l'equation de conservation de la masse multipliee par une 
composante de la vitesse a l'equation de conservation de cette composante de la quantite de mouvement) 



3v„ 



+ ^ = 



dx dy 



1 at 

dv 

It 



3v. 



+ pv x — + i n 

dv. 



dv^ = _dp 

y dy dx 

dv 



a 2 v„ a z v^ 



+ 



dx 2 dy 2 ) 



+ Pg* 



dx ^ y dy 



dp 
dy 



fd 2 V y d 2 V y ) 

d^^r^ 



Forme convective condensee 



Dv Dv 

V.v = p == -Vp + (lV 2 V + pg avec = p 

Dt Dt 


^ L - > 

+ — Vv 2 + ( V A V) A V 

v 3t 2 j 


= p 


— + (v.V)v 
v 3t J 


cartesien 
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Prise en compte des forces de gravite dans le cas d'un fluide incompressible 

Si les forces exterieures a distance derivent d'un potentiel, comme c'est le cas pour la gravite, le 
terme correspondant peut etre combine avec celui relatif aux forces de pression: 









^dz 1 dj} 




r 0) 


pV<£ = -pVgz 


= -Vpgz pour un fluide incompressible 


Vz = 


dz/dy 













\dzldzj 







d'ou -— + V.(pvv) = -V(p + pgz) + pvV 2 v 
at 

Dans le cas d'un fluide incompressible la pesanteur n'influe pas sur le champ des vitesses si le 
debit est impose. 

On note p =p+pgz la pression motrice (elle est parfois notee abusivement p) 
d'ou: + V.(pvv) = -Vp + pvV 2 v 

b) conditions aux limites du domaine 

Sur la section d'entree 1: la vitesse est connue: (v D , 0) 

surSl(x = 0, -Ll/2< y< Ll/2) : v x = v etv y =0 (eq.3) 

Sur la section de sortie 2: la pression motrice est connue, la vitesse est horizontale: 
sur S2 (x = L2 + L3, -Ll/2 < y < Ll/2) : p = p , v y = (eq.4) 

Sur les parois 3 a 10, le vecteur vitesse est nul: 

exemple: sur S3 (x = 0, LI/ 2 < y < L/ 2) : v x = 0, v y =0 (eq.5) 

c) Conditions initiales 

Nous allons considerez le cas ou le fluide est initialement au repos 
pour t=0 (quels que soient x et y): v x = v v .=0 (eq.6) 



Le probleme est ainsi entierement defini, nous avons combine 

- les equations de conservation (masse et quantite de mouvement) 

- les equations d'etat (p=cste) 

- les lois de transfert (fluide newtonien) 

- les conditions aux limites 

- et les conditions initiales dans un cas particulier. 
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2.6.2) Adimensionnalisation 



Definissons maintenant des variables sans dimension: 



x y 

- les composantes adimensionnelles du vecteur vitesse: v * = — , v * = — soit v* = v / v n 

v y v 

v * vaut au niveau des parois, et initialement dans tout le domaine 
v * vaut (1,0) sur la section d'entree 

- les coordonnees adimensionnelles: x* = x / L, y* = y / L 

En effectuant les changements de variables correspondant a cette adimensionnalisation on obtient: 
a) Interieur du domaine 



dw * 3v, 



- + - 



dx * dy * J 
3v„ * 



= (conservation de la masse, eq. 1 multipliee par L/vq) 

d 



P-Po 



v o- 1 

I L J 



+ 



3v *v * 9v *v * 



3x* 



dy : 



V P v o J 
dx * 



1 



+ 



V V 



'3 2 v x * a 2 v x *^ 

y dx* 2 dy* 2 J 



(conservation de la quantite de mouvement selon ox, eq.2 multipliee par — ; - ) 

PVo 

Par 1' analyse dimensionelle, on voit apparaitre un temps adimensionnel, ou L/v D est le temps 
necessaire pour parcourir la distance L a la vitesse vo . 



t* = 



v .t 



t 



L/v ^~ temps caractenstique 



On voit egalement apparaitre le rapport entre la vitesse a l'entree v Q et L/v que Ton appelle 
nombre de Reynolds (cf. nombre de Peclet pour les transferts de matiere) et que Ton note Re. 

On peut l'interpreter comme le rapport entre un flux convectif de quantite de mouvement (vitesse 
v ) et un flux diffusif de quantite de mouvement, c'est a dire une contrainte de cisaillement 
(gradient de vitesse v /L) 



Re = ^ 



p.v .v 



= 2 



pv(v /L) pv(v /L) 



pour un tube circulaire, on prendra le diametre pvD 

i Re = 

comme dimension caracteristique \i 



On peut encore considerer que le nombre de Reynolds compare la pression dynamique:l/2.p.v 2 
(energie cinetique par unite de volume) et une contrainte visqueuse de cisaillement. 

Enfin on voit que la pression peut etre adimensionnalisee par p.v 2 soit le double de la pression 
dynamique. La variable adimensionnelle correspondante est appelee nombre d'Euler. 



Eu = 



P-Po 
PVo 
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L'equation de continuite peut se reecrire: 



d\„ * 3v v * 



■ + ■ 



= soit vectoriellement: V* . v* = (eq. 1 ' ) 



dx* dy * 

ou V * represente le vecteur symbolique 



V 3x * dy *) 

L'equation de conservation de la quantite de mouvement peut se reecrire: 



d\ „ * 



■ + 



d\ *\ * dw x *v v * N 



at* 

vectoriellement 



dx* dy* 
d\ * 



3(Eu) J_ 
dx * Re 



a 2 v x * a z v x " 

dx * 2 dy* 2 



1 



+ V*.(v*v*) = -V*(Eu) + — V 2 *(v*) (eq. 2') 
3t* Re 



b) Limites du domaine 

LA 
2L 

L2 + L3 



surSl (x* = 0, - 
surS2 (x* = 



< y*< 



M ):v x * = l, v y * = 



— < y*<— ) : Eu = 0, v* = 
L L y 



(eq.3') 
(eq. 4') 



LI 



sur S3 (x* = 0, — < y* < 1) : v x * = 0, v * = (similaires pour S4 a S10) (eq. 5') 



c) Conditions initiales 

pour t* = 0, quels que soient x* et y*: v x * = 0, v y * = 



(eq. 6') 



Les equations 1' a 6' decrivent entierement le systeme. Elles permettent de calculer revolution du 
champ de la vitesse reduite: v * et de la pression reduite: Eu en fonction d'un seul parametre 
adimensionnel: le nombre de Reynolds et de la geometrie: (Ll/L), (L2/L) et (L3/L). 



(v*,Eu) = f{x*,y*,t*, Re, geometrie} 



Remarque : 



Pour un ecoulement incompressible (a debit impose), la pression n'intervient qu'au travers de son 
gradient. Le parametre p n'influe pas sur le champ des vitesses. Si Ton prend p ' a la place de p Q 
la pression sera juste decalee de (p' -p ) dans tout le domaine. 

Contrairement au cas du transport de matiere, le probleme n'est pas lineaire par rapport a la 
vitesse, car cette vitesse intervient au carre dans le terme de convection 

Si on forme alors la matrice des dimensions formee par les variables, on obtient. 





v x 


1* 


P 


X 


y 


t 


L 


v 


p 


v 


m 


1 


1 


-1 


1 


l 





1 


1 


-3 


2 


s 


-1 


-1 


-2 








1 





-1 





-1 


kg 







1 

















1 






Nous avons 10 parametres dimensionnels et le rang de la matrice des dimensions fondamentales 
vaut 3 (colonnes L, vq, p par exemples). La forme adimensionnelle de relation fait done intervenir 
(10-3)=7 parametres: (v x *,Vy*,Eu)=f{x*,y*,t*,Re) 
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2.6.3) Expression des equations de Navier- Stokes en coordonnees cylindriques 

Dans de nombreux cas pratiques on est amene a resoudre les equations de Navier stockes en 
coordonnees cylindriques. 

Le plus souvent (mais pas toujours) on peut considerer a priori une symetrie cylindrique: 
(3p / 30) = (3v r / 39) = (3v e / 30) = (3v z / 30) = 

Cela suppose que la geometrie, les conditions aux limites et les conditions initiales repondent a 
une symetrie axiale, mais la vitesse tangentielle ve peut ne pas etre nulle. Les equations se 
simplifient alors quelque peu: 

Hypotheses: symetrie axiale 
fluide newtonien 
incompressible 



3v r 3v„ 
■ + v. 



3t 



3r 

centrifuge 



+ V. 



3v^ 
3z 



J_3p 
p 3r 



+ v 



dfld, 0| 3 2 v 
K) + 



3r 



vr 3r 



3z 2 



3v a 3v e 



v.v c 



+ V. 



- + ■ 



3t r 3r _r_ 

Coriolis 



+ V. 



3Ve 

3z 



= V 



aria, ^ 3 2 v e 



3v, 3v 7 3v 



3t 



r 3r 



3z 



p 3z 



+ v 



1_3_ 
r 3r 



3v , l 3 2 v , 
+ ■ 



3r J dz 2 
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2.7) Notion de turbulence 



When I die and go to the heaven, there are two matters on which I hope for enlightenment. One is quantum 
electrodynamics and the other is turbulence in fluids. And about the former, I am really rather optimistic. 
Horace Lamb 

Reynolds a effectue une experience celebre, maintes fois renouvelee schematisee ci-dessous: 
^ colorant 




anemometre 



debit faible Re<2000 



1 



colorant 



anemometre 




debit eleve Re> 10000 

On opere avec des conditions aux limites quasiment permanentes, le debit est fixe a l'entree et 
la pression est maintenue constante a l'aval. 

Lorsque le nombre de Reynolds est faible, la vitesse est orientee suivant l'axe du tube, le profil des 
vitesses est parabolique comme l'indique le calcul (effectue en regime permanent), le colorant suit 
un filet fluide quasiment rectiligne comme on s'y attend puisque la vitesse radiale est nulle (il se 
disperse tres legerement du fait de la diffusivite moleculaire). Si on mesure la vitesse en un point 
avec un anemometre a reponse rapide 8 , on observe que la vitesse est constante dans le temps: 
l'ecoulement est dit laminaire.. 

Lorsque le nombre de Reynolds est eleve, la vitesse est en moyenne orientee selon l'axe du tube, 
mais le profil n'est plus parabolique (il est presque plat), le colorant diffuse tres rapidement comme 
si les trajectoires des particules de colorant etaient aleatoires et en tout cas beaucoup plus 
rapidement que ne laissent prevoir les equations de transport par diffusion moleculaire. Si on 
mesure la vitesse instantanee en un point, on constate qu'elle fluctue autour d'une moyenne, le 
regime d'ecoulement est done transitoire ! L'ecoulement est dit turbulent. 

Dans la pratique, pour une conduite circulaire, l'ecoulement est toujours laminaire pour 
Re=v d .D / v<2000, il est presque toujours turbulent pour Re>10000 (pour des installations 
industrielles on considere generalement que l'ecoulement est turbulent a partir de Re=4000). 
(on utilise dans la definition de Re, la vitesse debitante: vj qui est le debit divise par la section) 



Une telle mesure peut etre fait par exemple avec un anemometre a fil chaud. Un fil de platine de 5 microns de diametre est 
traverse par un courant electrique, il va done dissiper la chaleur produite par effet Joule dans le fluide. Pour une puissance 
donnee sa temperature sera d'autant plus elevee que la vitesse du fluide est faible. La resistivite electrique du platine varie 
avec la temperature, en mesurant la resistance on peut done connaitre la temperature du fil. 
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Dans le cas d'un ecoulement turbulent, on decompose la vitesse en une composante moyenne et 
une composante aleatoire: 

v x = v x + v x ' , v y = v y + v y ' soit vectoriellement; v = v + v' 

On peut alors reecrire les equations de Navier-Stokes et les moyenner dans le temps. Pour un 
ecoulement permanent en moyenne d'un fluide newtonien incompressible, on peut ecrire (quantite 
de mouvement dans la direction ox) 



3P(V X +V X ') 



at 



3p(v x +v x ')(v x +v x ') | 3p(v x +v x ')(v y +v y ') 



dx 



dy 



3(p + p') 



dx 



+ pv 



3 2 (v x +v x ') _ 3 2 (V x +v x ') 



dx : 



+ - 



dy 2 



On utilise alors les proprietes de linearite : 

- la moyenne de la derivee est egale a la derivee de la moyenne : 3(p + p' ) / 3x = 3(p + p' ) / 3x 

- la moyenne d'une somme est egale a la somme des moyennes : (p + p') = f +p' = p + 

On obtient ainsi une equation qui porte sur la vitesse moyenne qui ressemble beaucoup a 
1' equation en 1' absence de fluctuations de vitesse (Navier-Stokes en regime permanent laminaire) 
mais comportant un terme supplementaire qui provient du fait que : 

- la moyenne d'un produit n'est pas egale au produit des moyennes (v x + v x ')(v y + v y ') = v x .v y + v x '.v y ' 

C'est cette non-linearite qui explique l'effet papillon en meteorologie : un battement d'aile d'un 
papillon peut engendrer un cyclone (les fluctuations de vitesse ont une influence non negligeable 
sur le mouvement moyen). 



3 P V x V x | dp V x V „ f 



dx 



dy 



apv x 'v x ' + apVV 



dx dy 
terme supplementaire 



dp 

= -aI +pv 



r d 2 v, 3 2 v > 



dx 2 dy 2 



+ Pg, 



Le terme supplementaire fait intervenir la variance de v' x et la co variance de v' x et v' y . La matrice 
de covariance des fluctuations de vitesse multipliee par la masse volumique est appelee tenseur 
des contraintes de Reynolds . Les composantes de cette matrice ont la dimension d'une densite de 
flux de quantite de mouvement c'est a dire d'une contrainte (Pa). Ces contraintes sont dues a 
1' agitation turbulente (a l'echelle des tourbillons), de meme que les contraintes visqueuses sont 
dues a l'agitation thermique (a l'echelle moleculaire): 



Tenseur des containtes de Reynolds: pv'v' = p 



v v v v 

X x X y 



\ V V V V / 

V y x y y / 



Les equations du mouvement moyen s'ecrivent alors sous forme condensee 



V.V = et V.(pW) = -Vp + pvV 2 V-V.(pv'v') + pg 
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Pour traduire les flux de quantite de mouvement dus au mouvement turbulent, on utilise le plus 
souvent la notion de viscosite turbulente. Ainsi par analogie on ecrira : 



^xy = -PV 



— - + — - 
dy dx 



= -2pv 



dx 



x = -pv(Vv + (Vv) t ) 

viscosite cinematique 

v : 

moleculaire (ou laminaire) 



ana log ie 



ana log ie 



-> Pv x 'v y ' = -pv t 



3v„ dv x 



+ ■ 



* Pv x 'v x ' = -2pv t 



3y 3x 



dv. 



3x 



p v'v' =-pv t (Vv + (Vv) t ) 

viscosite cinematique 



turbulente 



Contrainte visqueuse : agitation moleculaire 

A 



Contrainte turbulente : fluctuation de vitesse 



En adoptant un modele de viscosite turbulente, tout se passe comme si on pouvait utiliser, en 
ecoulement turbulent, les equations de Navier Stokes en regime permanent en rempla§ant la 
viscosite par la somme de la viscosite moleculaire et de la viscosite turbulente. 

V : laminaire — > (v + v t ) : turbulent 



V.(pvv) = -Vp + pV. 



(v + v t )(w + (VV)') 



+ Pg 



V.v = 



II faut remarquer que dans ce cas V t depend de l'ecoulement local. Ce n'est pas uniquement une 

caracteristique du fluide comme v. Toute la difficulte consiste a modeliser alors la viscosite 
turbulente, ceci sort du cadre de ce cours. 



La forme adimensionnelle des equations de Reynolds s'ecrit : 

Re. V* .(V * V*) = -Re. VEu + V 2 V * - Re. V.(v'*v'*) V* . V* = 

On voit que l'importance relative de la turbulence est d'autant plus grande que le nombre de 
Reynolds est eleve. C'est l'inertie du fluide qui provoque la turbulence, par contre les frottements 
visqueux l'amortissent. Le nombre de Reynolds compare l'ordre de grandeur de l'energie 
cinetique volumique : l/2pv 2 et l'ordre de grandeur des contraintes visqueuses: pv.Vq/L, il est 
done un bon indicateur du type d'ecoulement (laminaire ou turbulent). 
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2.8) Equations du mouvement pour un fluide parfait 9 



Un fluide est dit parfait si les forces visqueuses peuvent etre considerees comme negligeables 
devant les forces d'inertie 10 . Dans ce cas, les contraintes totales se reduisent aux contraintes de 
pression. L'expression de la conservation de la quantite de mouvement sur un element de masse 
fixe suivi dans son mouvement s'ecrit alors 



1 i 



Dv 
P ^t =P 



^ ||v(v.v) + (Vav)av 



.at 2 

Masse volumique constante / Forces a distance = gravite 

Dv 



= -V.7T + pg = -Vp + pg 



P ^T = P 



av i . /- \ 

— + -V(v.Vj+ V A V A V 

at 2 v ; 



-V(p + pgz) 



av 
at 



+ (V A v) A V 



= - V (P+2p V + P§ Z ) 



Projetons cette equation sur une ligne de courant pour un ecoulement permanent (la projection se 
traduit par un produit scalaire avec le vecteur vitesse) 



3v 

at 



= et (Vav)av.v = => V 



p + -pv 2 + pgz|.v = 



Pour un ecoulement permanent d un fluide parfait, la charge reste constante le long d'une 
ligne de courant. 

Si l'ecoulement est permanent et irrotationnel (Vav = 0), la charge est constante dans tout 
l'ecoulement . 

Si l'ecoulement est irrotationnel, on peut demontrer que le champ des vitesses derive d'un 
potentiel: v = . La conservation de la masse s'ecrit alors V.v = V^V 1 !') = V 2 ^ = . II suffit 
dans ce cas, de trouver une fonction *P de Laplacien nul qui verifie les conditions aux limites. 

Masse volumique fonction de la pression / Forces a distances derivant d'un potentiel 4> 

La relation precedente se generalise sous la forme : 



3v 

at 



+ (V A v) A V = -V| 



i 



dp' 
P op{p'} 



+ — + $ 

2 



par exemple 4> = gz 



9 Rappel : un fluide parfait etant sans frottement, sa vitesse tangentielle le long d'une paroi n'est, a priori, pas nulle (la 
condition aux limite d'adherence n'est pas remplie). Par contre la vitesse normale a la paroi est nulle. 

10 On devrait plutot parler d'ecoulement parfait de fluide. L'hypothese negligeant les forces visqueuses est plus ou 
moins verifiee pour tel ou tel type d'ecoulement (pour des Reynolds eleves). Selon le cas, l'eau ou l'air peuvent etre 
consideres comme fluide parfait ou fluide visqueux. 

11 En developpant on obtient les equations d'Euler (voir equations de Navier Stokes avec u=0) 
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2.9) Bilan de quantite de mouvement sur un systeme macroscopique ouvert 



Pour effectuer le bilan des forces sur le fluide qui traverse un systeme ouvert, il suffit d'effectuer 
une integration sur le volume considere. Pour cela on utilise, comme au paragraphe 1.7, la formule 
de Leibnitz et le theoreme d'Ostrogradskii. Pour la resultante, on obtient 12 : 



d 

— j pv x dV 
dt J v r x 


+ | s Pv x (v-v s ). ndS = 


-J s Pn.e x .dS 


- f x „ £. dS 


+ J v p gx dv 


d f 

— pvdV 
dt J v^ 


+ J s p(v(v-v s )). ndS = 


-J s pndS - 


f x. ii dS + 


J v pgdv 


accumulation 


bilan des flux convectifs 


contraintes 


contraintes 


forces exterieures 






de pression 


visqueuses 

V 


a distance 






resultantes 


des forces 


exterieures 



Le terme p(v - v s ). n dS represente le flux de matiere (masse totale) sortant a travers une section 
elementaire dS de vecteur normal unitaire exterieur ii (cf. § 1.7.1) . Le terme pv(v-v s ). n dS 
represente le flux de quantite de mouvement transporte par ce flux de matiere (flux convectif). 
Cette expression est similaire a celle du flux convectif du constituant i obtenue au § 1.7.2: 
p©i(v-v s ). ndS. 

Les termes ( - p ii dS ) et ( - x. ii dS ) representent la force de pression et la force visqueuse 
exercee sur le fluide contenu dans le systeme a travers une section elementaire n dS 

Si Ton s'interesse a la composante selon ox de la force visqueuse a travers une paroi inclinee 
(figure de droite ci-dessous) on peut faire un bilan sur les faces du triangle rectangle 
correspondant. La face horizontale (_L a oy) a une surface egale a (-n y dS ). Pour avoir la force dans 
la direction ox, on multiplie cette surface par la contrainte x yx . On fait de meme pour la surface 
verticale et on prend la somme, on obtient alors : -(x xx .n x +x yx .n y ) dS , c'est a dire la 

composante selon x de - x. ii dS . 




-n y dS x. 



Le systeme peut avoir des frontieres mobiles se deplacant a la vitesse v s 
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On peut ecrire egalement une relation equivalente pour le moment par rapport a un point 0. Ce 
point peut ne pas appartenir au domaine. Pour ecrire le moment par rapport a un axe passant par 0, 
il suffit de projeter la relation sur l'axe, c'est a dire la multiplier scalairement par le vecteur unitaire 
de l'axe. 



j t I p(x a v)dV + J s p((x a v)(v - v s )). n dS = 


- J s p(x a n)dS 


- j" s XA(T.n)dS + J v p(xAg)dV 


accumulation bilan des flux convectifs 


contraintes 


contraintes forces exterieures 




de pression 


visqueuses a distance 




moment 


des forces exterieures 



Pour un systeme ayant une section d' entree et une section de sortie et en negligeant les contraintes 
visqueuses sur ces sections on obtient ainsi le theoreme d'Euler : 



dP / - _-\ 

— + A mv + pS) = mg+ F e; 

at s-e 




(pS) e 

(mV) e 



P est la quantite de mouvement du fluide contenu dans le domaine: P = | pvdV 



v est la vitesse de melange: v = 



fpv((v-v s ).n)dS 



m 



si la section d' entree est fixe, si la direction de la vitesse est identique sur toute la section 
d' entree et que p est uniforme sur cette section, alors v e est le vecteur oriente selon cette 

direction et de norme egale a J s vys / J s vdS 

si la vitesse est uniforme, on a simplement: v = v 



S e est le vecteur de norme egale a la section d' entree, normal a cette section et oriente dans 
le sens du flux. 

p e est la moyenne surfacique de p sur la section d' entree (qui est en general uniforme) 

F ext est la resultante des forces (sauf celles dues a la gravite et a la pression sur les sections 
d' entree et de sortie), il s'agit typiquement des forces exercees par les parois fixes ou 
mobiles (pales d' une pompe par exemple) du systeme sur le fluide qu'il contient. 
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Pour un systeme ayant une section d'entree et une section de sortie, en negligeant les contraintes 
visqueuses (flux diffusifs) sur les sections d'entree et de sortie et si la vitesse, la pression et la 
masse volumique sont uniformes sur chacune des sections d'entree et de sortie: 



dn 







A (rh(x a v) + p(x a S)) = m(0Mg a g) + M 



+ 



O.ext 



dt 



n o est le moment par rapport a 0, de la quantite de mouvement du fluide contenu dans le domaine: 
n o = I p(x a v)dV (x est le vecteur qui relie au point considere du domaine) 



v e est la vitesse (supposee uniforme) sur la section d'entree 

x e est le vecteur qui relie au centre de la section d'entree 

(on suppose que |x e | est grand devant les dimensions de la section d'entree) 

M g est le centre de gravite du fluide contenu dans le systeme 

M ext est le moment par rapport a des forces exterieures (sauf celles dues a la gravite 
et a la pression sur les sections d'entree et de sortie) 

Pour un systeme tournant autour d'un axe (passant par O), on projetera cette equation sur le vecteur 
unitaire de cet axe (si l'axe est vertical, cela revient a prendre la composante selon z de 1' equation) 

2.10) Bilan d'energie mecanique sur un element de volume fixe 

En mecanique du point materiel le travail (resp. la puissance) d'une force est egale au produit 
(scalaire) de la force et du deplacement (resp. de la vitesse). 

Exemple: 

Si un wagonnet monte une cote a vitesse constante ( d(mv)/dt= ) alors les forces exterieures 
(reaction des rails: Fr, forces de poussee: Fp et forces de gravite: Fg) s'equilibrent. La puissance de 
la force de poussee est egale a: 



W p =F p .v 
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De la meme facon, pour obtenir les puissances des differentes forces exercees sur un fluide, nous 
multiplions scalairement l'equation de quantite de mouvement par le vecteur vitesse. Nous 
supposerons que les seules forces a distance sont celles de gravite. 



v.^^ + v.(v.(pvv)) = pv.g- v.Vp- v.(v.t) 



et g = -Vgz 



3pv x _ 9pVj 



3t" +v ^ at 



3(pv x v x ) 3(pv y v x ) 



^3(pv y v x ) 3(pv y v y )^ 



dx 



dy 



+ v. 



dx 



dy 



J J 



+ p(g x v x+gy v y )-(v x ^ + v y ^ 



V 

V v 



dx dy 



3(pT yx ) + 3(px yy ) 



dx dy 



J J 



Apres un calcul fastidieux mais qui peut etre traite de fa§on elementaire on obtient une equation de 
bilan sur un element de volume fixe dans l'espace de l'energie mecanique (cinetique et potentielle 
de gravite) 13 . 



3(ipv 2 ) 

at 



+ v 



.(^pv 2 vj accumulation et bilan (s - e) des flux convectifs d'energie cinetique 



3(pgz) 

+ — - — +V.(pgzv) accumulation et bilan (s-e) des flux convectifs d'energie potentielle 
at 

= (- V.(pv)) puissance fournie (par unite de volume) par les forces de pression 

(pour un fluide incompressible, la pression peut etre consideree comme 
une energie par unite de volume, ce terme represente alors le bilan (e - s) 
des flux convectifs "d'energie de pression" ) 

- (- p(V.v)) flux d'energie mecanique convertie de fagon reversible en energie interne 

(par unite de volume) (terme nul pour les fluides incompressibles) 

+ (- V.(x.v)) puissance fournie (par unite de volume) par les forces visqueuses 

- (- x:(Vv) ) flux d'energie mecanique convertie de facon irreversible en energie interne 

(par unite de volume) : dissipation sous forme de chaleur 

(d' apres le second principe de la thermodynamique ce terme est toujours positif) 



Cette equation n'est pas une equation fondamentale mais a ete derivee des precedentes. 



L'equation peut etre generalised dans la mesure ou les forces volumiques derivent d'un potentiel constant dans le 
temps (forces electrostatiques, forces centrifuges, etc. en plus de la force de gravite) 
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2.11) Bilan d'energie mecanique sur un domaine macroscopique 



Nous allons maintenant ecrire le bilan d'energie mecanique pour un fluide traversant un systeme 
macroscopique. 

- en regime permanent: 3/3t=0 

- pour un fluide incompressible: p=cst => V.v = 

- lorsque les forces sont celles de gravite (elles derivent d'un potentiel) 14 : g = -Vgz 

- en negligeant la puissance des forces visqueuses sur les sections d'entree et de sortie 

- sans apport d'energie mecanique: W = 

(il n'y a pas de pieces en mouvement dans le fluide comme ce serait le cas dans une pompe par ex.) 

En integrant l'equation du paragraphe precedent sur le volume de fluide contenu dans le systeme, on 
a: 

J V.(p(^ v 2 + gz)v)dV = -J V.(pv )dV - J V.(f . v)dV - J (-f I Vv)dV 

V V V V 

En utilisant le theoreme de Gauss Ostrogradskii, on obtient 

J p(^ v 2 + gz)(v . n)dS = - J p(v . n)dS - J (f . v) . ndS - J (-^ Vv)dV 

S S S V 

La vitesse est nulle sur les parois et on suppose que (x.v) est negligeable sur les sections 
d'entree - sortie (assez bien verifie en general), d'ou: 

J(^pv 2 +pgz + p)(v.n)dS = - J(-^:Vv)dV 

Se+Ss V 

Si par ailleurs on apporte une puissance mecanique W ext au fluide contenu dans le systeme (par 
l'intermediaire de parois mobiles : pales d'une pompe par exemple), alors: 



J(ipv 2 +pgz + p)(v.n)dS 


- J(-x:Vv)dV 

V 


+ w ext 


Se+Ss 




bilan des flux d' energie 


energie mecanique dissipee 


energie mecanique 


cinetique, potentielle 


en chaleur du fait des 


effectivement apportee 


et de pression 


frottements visqueux 


au fluide contenu dans 


( sortie - entree ) 


(terme toujours negatif) 


le systeme 



regime permanent, fluide incompressible 



II suffit en fait que les forces a distance derivent d'un potentiel constant dans le temps ( 4><->gz) 
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On peut ecrire de facon simple le bilan d' energie mecanique 



pour un systeme a une entree et une sortie 

en regime permanent 

pour un fluide incompressible 

■ lorsque les forces volumiques se 
resument a la gravite 
si sur la section d'entree et de sortie la 
pression motrice p+pgz est uniforme 



(Vp*) e 




(Vp*) s 



A (v.p *) = W ext - W d avec p* = p + ^ p v 2 + pgz : charge = energie mecanique par unite de volume 



W ext est l'energie mecanique apportee au fluide contenu dans le systeme par unite de temps : 

( puissance des forces exterieures sauf celles de gravite et de pression sur les entree-sortie, 
typiquement il s'agit de la puissance des forces exercees par les pales d'une pompe par ex. ) 

W d est l'energie mecanique dissipee en chaleur par unite de temps du fait des frottements 
visqueux dans le fluide contenu dans le systeme. 

^v 2 est l'energie cinetique (par unite de masse) de melange: \ v 2 = j\^pv 2 (v.n)ds/m 

- si sur chacune des sections d'entree et de sortie la vitesse possede une direction 
uniforme, alors: v 2 = j s v 3 dS j dS 

- si la vitesse est uniforme, on a simplement: v 2 = v 2 

Nous exprimons ici la charge : p* = p + ^pv 2 + pgz en Pascal, mais certains auteurs la divisent 

par pg et obtiennent ainsi une expression homogene a une hauteur. Elle represente la somme des 
energies cinetiques, potentielles et de pression 15 par unite de volume. 

Si Ton peut negliger les frottements visqueux (hypothese des fluides parfaits), la charge reste done 
constante au travers d'un systeme ou le long d'un tube de courant (hyp. 3/3t=0, p=cst, W ext =0) 

L'expression precedente peut etre generalisee au cas ou la masse volumique n'est pas constante 
mais depend seulement de la pression, et pour d'autres energies potentielles que celle de gravite, 
on raisonne alors par unite de masse: 



dp W ext 



W, 



m 



d 

m 



(regime permanent) exemple : 4> = gz 



Pour un fluide newtonien, on peut expliciter l'energie mecanique dissipee en chaleur par 
l'intermediaire des forces visqueuses. Peu importe pour nous cette expression, retenons 
simplement qu'il est possible de la calculer. 



W d = J v ( _ T - ^) dV = P V lv ((^ + ' ) : ^) dV = pV lv *^ dV 



avec ^> . 



dx 



+ 



Uy 



+ 



^3v_ 



dz 



+ 



+ 



3x dy 



+ 



3v. 3v x 



+ - 



dy dz 



+ 



3v„ 3v, 
+ 



dz dx 



15 Ce n'est que pour un fluide incompressible que Ton peut assimiler la pression a une energie par unite de volume 
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En regime permanent et pour une masse volumique constante (hypotheses enoncees 
precedemment) on a:V e =V s =V. On peut done diviser l'expression du bilan d'energie 
mecanique par le debit volumique: 

P*-p!=(W ext /v)-(w d /v) en J.m 3 = Pa 

p s - p e est la variation de l'energie mecanique par unite de volume du fluide entre l'entree et la 
sortie du systeme 

W ext /V represente l'energie mecanique apportee effectivement 16 au fluide par unite de volume au 

niveau d'une pompe ou d'un ventilateur par exemple 
W d /V represente la dissipation d'energie mecanique par unite de volume traversant le systeme 

du fait des frottements dans les tuyauteries d'un circuit par exemple, ce terme est appelee 
perte de charge 

Si le systeme considere est un circuit ouvert comprenant des tuyauteries, des coudes, des vannes 
etc. et eventuellement une pompe ou un ventilateur, on peut ecrire : 

Ps -Pe = A P pompe " A P circuit I 




p s - p e est la variation de l'energie mecanique (par unite de volume) du fluide entre l'entree et la 
sortie du systeme 

A P pompe est l'energie mecanique nette apportee par la pompe au fluide par unite de volume, 

ce terme est appele elevation manometrique de la pompe 
A P circuit est l'energie mecanique dissipee en chaleur du fait des frottements dans les differents 

elements du circuit (par unite de volume de fluide traversant le circuit), 

ce terme est appele perte de charge du circuit 

L'equation de bilan signifie simplement que la variation d'energie mecanique entre l'entree et la 
sortie du systeme est egale a l'energie mecanique apportee diminuee de l'energie mecanique 
dissipee par frottement 



Une partie de l'energie apportee au niveau de Farbre d'une pompe peut etre dissipe en chaleur dans les roulements 
et les joints 
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Dans l'application de tout bilan, c'est a vous de definir les limites du systeme, sachant que vous 
aurez a calculer ou exprimer certaines grandeurs aux frontieres du systeme. Pour un bac de 
lancement (resp. de reception), il est souvent commode de choisir comme entree du systeme (resp. 
comme sortie) un plan horizontal situe juste en-dessous du niveau du liquide car il est facile d'y 
exprimer la charge (cf. schema ci-dessus). Une alternative consiste a considerer une frontiere 
situee juste avant le retrecissement (resp. juste apres l'elargissement) ou l'energie cinetique du 
fluide peut encore etre consideree comme negligeable (cf. schema ci-dessous). 



a 




Remarque: Ap pompe est la variation de la charge (energie mecanique par unite de volume) entre 

l'entree et la sortie de la pompe, il s'agit de l'energie mecanique transmise au fluide par les pales 
diminuee de la dissipation d'energie mecanique par frottement visqueux a l'interieur de la pompe. 
La puissance utile de la pompe (puissance mecanique nette transmise au fluide) est simplement 

egale a V.Ap* mpe 

Le bilan d'energie mecanique sur un circuit ouvert peut encore s'ecrire : 



Ap 



pompe 

caracteristique 
pompe 



Ps -Pe + A Pcir, 

caracteristique 
reseau 



Ceci peut se representer graphiquement en faisant apparaitre le point de fonctionnement comme 
1' intersection de la caracteristique de la pompe et de la caracteristique du reseau qui dependent 
tous deux du debit: 



<- 



caracteristique reseau : p s - p e + Ap cin 
point de fonctionnement 



caracteristique de la pompe : Ap 



pompe 



Le cas des circuits fermes peut etre considere comme un cas particulier ou l'entree et la sortie 
sont confondues en un point du circuit, dans ce cas on a done 

pompe 

= Ap 

circuit 

Ceci signifie que l'energie mecanique apportee au niveau de la pompe compense la dissipation par 
frottement dans le circuit. 
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2.12 ) Caracteristique des pompes et ventilateurs 



La caracteristique essentielle d'une pompe est la relation qui lie (pour un fluide et une temperature 
donnee: v,p fixes ) l'augmentation de la charge entre l'entree et la sortie (appelee elevation 
manometrique) au debit de fluide; ceci en fonction de la vitesse de rotation. On peut ecrire: 



A p* = f{v,N,p,v,D} A!! = fj^.,m 

pompe pN 2 D 2 [ND 3 U. 



A p* Pompe type xxx 
eau 20 C 




2.13) Calcul des pertes de charge 

Pour calculer les pressions, les debits ou les puissances mecaniques pour un systeme donne, il 
reste a pouvoir calculer les pertes de charges notee simplement Ap* dans ce paragraphe. On 
distingue deux cas: 

Pertes de charge lineaires 
Le fluide s'ecoule dans une direction donnee, a travers une section constante et avec un profil 
de vitesse constant (etabli). Exemple: troncon de tuyauterie. Dans ce cas la perte de charge 
varie proportionnellement (de facon lineaire) avec la longueur de l'element: L. 
- Pertes de charge singulieres 

Exemple: coude, elargissement, retrecissement, vannes... 

a) Pertes de charge lineaires 

La perte de charge par unite de longueur depend de 5 variables (et de la geometrie du canal) 

re— L ^ 




on note pour simplifier v la vitesse debitante: v = V / S 
£ est la rugosite du tuyau (hauteur moyenne des asperites) 
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Trois unites fondamentales sont en jeu (kg,m,s), l'analyse dimensionnelle montre alors que la 
relation precedente (reliant 6 variables dimensionnelles) peut etre simplified en une relation entre 3 
nombres sans dimension. (Ap*/L est considere comme une seule variable) 

fcoeff. de. perte. de^j f nombre rugositel . _ Ap * fvD el 

V charge, lineaire J ~ [de. Reynolds relative J ^ f_ Ll 9 _ 1 v ' D I 

pv L J 

D 2 



Le nombre de Reynolds peut etre exprime en fonction 

- de la viscosite cinematique: v : Re 

- ou de la viscosite dynamique: (l=pv : Re 

Attention: On utilise egalement le nombre de Fanning : Fa : 
Conduites circulaires 

64 

Pour une conduite circulaire, en regime laminaire (Re<2000) on a: Af = ^ 
En regime turbulent (Re>4000), differentes correlations sont utilisees 

- citons celle de Colebrook pour les tubes rugueux qui est representee par l'abaque de Moody. 



vD 
v 

pvD 

^ rz.paroi ^ f 



1 

= -2 log! 



251 

■ + ■ 



^ 5l \3.71D Rejk^, 

- et celles Blasius pour les tubes lisses 

A f = 0.316 Re 1 ' 4 pour 4000 < Re < 10 5 et X f = 0.18 Re" 1 ' 5 pourRe>10 5 
Conduites non circulaires 

En regime turbulent on peut utiliser les memes relations pour des conduites non circulaires en 
rempla§ant le diametre par le diametre hydraulique (pour une canalisation partiellement remplie 
d'eau, on prendra la section de passage du liquide et le perimetre mouille c'est a dire la longueur 
ou s'exerce le frottement liquide/paroi, on neglige le frottement liquide/air ) 

4. Section 
^ Perimetre 

Attention: Certains auteurs utilisent egalement le rayon hydraulique, mais celui-ci 
ne correspond pas au rayon pour une conduite circulaire: a eviter ! 

b) Pertes de charge singulieres 

C'est ici la perte de charge (et non plus la perte de charge par unite de longueur) qui depend de 5 
variables (et de la geometrie du systeme). 

Ap* = f{ v, D, 8, p,v} 
L'analyse dimensionnelle aboutit alors a une relation entre 3 nombres sans dimension: 



j^coeff. de. perte. de^j _ Ap* fvD e 
Uharge.singuliereJ _ ^p v 2 _ [ v 'D 
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Dans le cas de coudes, vannes etc., on definit en general une longueur equivalente de tuyauterie 
circulaire: L e . Par definition, une tuyauterie de cette longueur (et dont le diametre correspond a 
celui de l'obstacle: vanne...) genererait la meme perte de charge que l'obstacle. On ecrit alors: 

Ap* _ JvD 8 1 

D 2 V 

Le terme de droite est le coefficient de perte de charge lineaire de la tuyauterie cylindrique raccordee 
a la vanne, au coude etc. La valeur de L e /D ne depend que des rapports geometriques et non pas de 
la dimension. Par exemple pour un coude a 90°, L e /D est fonction du rapport entre le rayon de 
courbure du coude et le diametre du tube. 



Dans le cas de retrecissements ou d'elargissements, on neglige en general les pertes par frottements 
visqueux ainsi que les effets de rugosite, et on ecrit: 



ou K est une constante qui ne depend que de la forme du retrecissement 
(elargissement) et du rapport des sections d'entree et de sortie 

par convention on prend la vitesse dans la section la plus petite 



pour un elargissement au niveau de l'embouchure d'une conduite dans un grand bac K=l 

L'utilisation de longueurs equivalentes reduites (L e /D) et de coefficients de perte de charge 
singuliere independants du Reynolds n'est valable qu'en premiere approximation. Pour plus de 
details et pour des applications numeriques on se reportera a des ouvrages specialises tel que le 
Memento des pertes de charges (Idel'cik I.E. Eyrolles 1979), on trouvera egalement quelques 
indications dans Mecanique et rheologie desfluides en genie chimique (Midoux N. Lavoisier 1988). 



2.14) Calcul de circuit 



Dans le cas d'un circuit compose de differents sous systemes (tubes, vannes, reservoir) et 
remplissant les hypotheses suivantes: 

-regime permanent 

-fluide de masse volumique constante et newtonien 

-pression (et potentiel) uniforme sur les sections de raccordement des sous-systemes 
-pour chaque section tubulaire la longueur d'etablissement est negligeable devant la 
longueur du troncon (distance au bout de laquelle le profil de vitesse est constant) 
-un seul diametre de tubulure 
-une seule rugosite de paroi 



* * 1 2 

Ap* = -pv 2 

circuit L 




y D D y 




avec. X f = f{Re,8/D} 



Si l'une des hypotheses ci-dessus n'est pas verifiee, remonter en sens inverse les differentes 
etapes jusqu'a F=my si necessaire ! 
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formule valable 



formule non valable ( longueur d'etablissement non 

negligeable sur le troncon vertical) 



(la formule donnerait d'ailleurs le meme resultat dans ces 2 cas 
alors que la perte de charge est plus forte a gauche qu'a droite) 

Cette relation (eventuellement adaptee si la situation plus complexe que celle decrite pas les 
hypotheses d' application) est a combiner dans la pratique avec l'equation de bilan d'energie 
mecanique deja enonce precedemment 

pour un circuit ouvert (en regime permanent pour un fluide incompressible) : 

Ps -Pe =Ap pompe -A P 

circuit 

pour un circuit ferme (en regime permanent pour un fluide incompressible) : 



Ap 



pompe 



2.15) Ecoulement autour d' obstacles 

Considerons un objet (obstacle) place dans un ecoulement de direction principale Ox. La 
composante selon Ox de la resultante des forces, appelee force de trainee, est fonction de la 
geometrie ainsi que : 

- d'une dimension caracteristique: D 

- de la masse volumique du fluide: p 

- de la viscosite cinematique du fluide: v 

- de la vitesse du fluide (orientee selon Ox) en amont de l'objet: v 
F x =f{D,p,v,v} 

Cette relation fait intervenir 5 variables qui s'expriment a partir de 3 unites fondamentales: (m,s,kg) 
(le rang de la matrice des dimensions est egalement 3). L'analyse dimensionnelle nous indique alors 
que Ton peut exprimer cette relation a l'aide de (5-3)=2 nombres sans dimension: 

C x = f {Re} avec C x = et Re = — 

yP v2 v 

ou A est la section droite de l'objet perpendiculairement a la direction de l'ecoulement 
(pour une sphere A=7lD 2 /4, ne pas confondre avec la surface de l'objet: S=7tD 2 ) 



Cas d'une sphere: Re<2 C x =24/Re Stokes 



2<Re<500 C x =10/VRe Lapple et Sherpherd 

(il existe des relations plus precises) 

Re>500 C x =0,43 Newton 



Si l'objet n'est pas symetrique par rapport a l'axe Ox, il peut egalement y avoir une composante selon 
Oy et Oz (portance) et un moment des forces non nul par rapport au centre de l'objet. 
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2.16) Ecoulement en milieu poreux: perte de charge a travers un lit fixe 

On considere que l'ecoulement se fait a travers un canal de section tres complexe mais qui peut etre 
caracterise par son diametre hydraulique moyen . 

sec tion occupee par le fluide 

D, =4 — 

perimetre mouille 

Integrons numerateur et denominateur entre et L 

volume occupe par le fluide 
h surface des particules 

Divisons par le volume du lit et definissons £ comme le taux de vide, on obtient alors: 

4.8 surface des particules 

D. = oil a v est la surface specifique des particules: 

a v (l-8) ' volume des particules 

Pour des spheres, on a: a v = (tcD 2 ) / (tiD 3 / 6) = 6 / D 

Pour des particules non spheriques on definit un diametre equivalent D = 6 / a v 

Ap * k 

En regime laminaire, on sait que: — = — (k = 64 pour des tubes circulaires) 

pv 2 ^ e 

Ap* u.v (1-e) 2 

Si on applique cette relation aux lits fixes, on obtient: = 72 — 

L D 2 e 3 

Mais la tortuosite des canaux fait que le coefficient 72 est sous-estime, on trouve en fait 
(experimentalement) : 

Ap* = 15() u..v(l-e) 2 



3 



D 2 8 

En regime turbulent, lorsque le nombre de Reynolds tend vers l'infini, la viscosite n'intervient plus 

Ap* 



L 1 2 

pv 2 

D h 2 K 



tend alors vers une constante 



Pour de grands nombres de Reynolds, on trouve experimentalement: 

Ap* = 1 75 pv 2 (l-e) 

L D 8 3 

La relation d'Ergun combine les deux relations precedentes: 



^! = 150 t^ilzi)! +1 , 75 P v2 (i-« 



D 2 e 3 D e 3 



v est la vitesse debitante: debit du fluide/ section totale du lit 

D est le diametre des particules (equivalent si les particules ne sont pas spheriques) 

8 est le taux de vide 

Le premier terme est preponderant pour de faibles nombres de Reynolds de particule : (vD)/(ev), 
le second terme predomine pour de grands nombres de Reynolds. 
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3) Transferts d'energie 



// vaut mieux acheter un thermometre en ete qu 'en hiver. II coute le meme prix mais il contient 
davantage de mercure. Francis Blanche 

Nous allons analyser maintenant le transport d'energie et plus particulierement de la chaleur dans 
un fluide dont le mouvement est connu. Dans un premier temps nous considerons un petit element 
de volume parallelepipedique (Ax.Ay.Az) fixe. 

3.1) Equation de conservation sur un element de volume fixe (energie totale, energie interne) 
3.1.1) Energie totale 

Le premier principe de la thermodynamique indique que 1' accumulation plus le bilan des flux 
d'energie interne et d'energie cinetique est egal a la somme du flux de chaleur et de la puissance de 
toutes les forces exterieures. 



d(U + K) 
dt 



+ A rh(U + K) = Q + W tc 



Ecrivons maintenant ce bilan sur un petit element de volume, nous aurons 

3p(u+» 



+ V.(p(u + »v) = -V.(jJ 



+ r„ 



+ pg.v 



bilan des 



dt vrv 2 ' I yJq > q 

puissances des 

accumulation bilan des flux bilan des flux source volumique puissance des forces exercees a 
convectifs diffusifs de chaleur de chaleur forces de gravite travers les facettes 

(effet Joule...) 

La puissance d'une force s'obtient en la multipliant scalairement par la vitesse, c'est ce que nous 
avons fait pour la gravite. Calculons maintenant le bilan des puissances des forces exercees a 
travers les facettes. 

Ax 



■xy 



(x-Ax/2 , y) 




La force exercee par le fluide situe a gauche sur celui de l'element vaut : Ay.Az.7C x> = Ay.Az. 
Cette expression devant etre prise en (x-Ax/2,y). 

La puissance de cette force vaut : Ay.Az.S x> .v = Ay.Az. [n xx .v x + 7C xy .v y ) toujours en (x-Ax/2,y) 
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Pour la puissance de la force exercee par le fluide situe a droite sur celui de 1' element, il convient 
de mettre un signe negatif (loi de Taction et de la reaction) car par convention la contrainte 7t x> 
est celle exercee dans le sens des x croissants (par la gauche sur la droite). Le centre de la facette 
est cette fois situe en (x+Ax/2,y). 

Le bilan sur les facettes de gauche et de droite divise par le volume de l'element s'ecrit done : 
(^xx-v x +^ xy .v y ){x-Ax/2,y}-(7C xx .v x +7t xy .v y ){x + Ax/2,y} ^( Ic n ' v ,+ Jt iy v yj 



Ax 



->-- 



On trouve de la meme facon pour le bilan sur les facettes du haut et du bas 
Finalement on peut condenser l'ecriture du bilan sous la forme suivante : 

5 Kx- V x +7 V V y) a Kx- V * +7 V V y) 



M ^ D "' dx 

3 (?V v x +7 V v y) 



3y 



dx 



■ + ■ 



= -V.(ji.v) 



En effet Tt.v = 



^xx ^xy 
V^yx ^yyj 



^xx- V x +^xy V y^ 



produit d' une matrice et d' un vecteur 



D'OU V.(tT.v) = 



fd/dx" 




la/ayy 


• 

V 



7t„.v„ + 7t„„.v 
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yx 



V. + 7t,„,.V y y 



xy y 

yy 



V + 7t .V 

x ' " xy " y 



V +71 .V 

x '"yy y, 



dx dy 
produit scalaire de deux vecteurs 



Ainsi le premier principe de la thermodynamique ecrit sur un petit element de volume (systeme 
ouvert) s'ecrit : 





v) = -V.(j q ) 


+ r q 


- V. (jt.v) 


+ pg»v 


accumulation bilan des flux 
convectifs 


bilan des flux 
diffusifs de chaleur 


source volumique 

de chaleur 
(effet Joule...) 


bilan des puissances 
des forces exercees a 
travers les facettes 


puissance des 
forces volumiques 
a distance 



On rappelle que la divergence V.(|) represente le bilan (sortie - entree) de la quantite dont la 
densite de flux est representee par le vecteur <j) (divise par le volume de l'element). 



Comme pour un fluide on a: n = pi + x , on peut ecrire: 



3 ^^ + V. (pfu + lv^v) =- V.(j q ) + r q - V. (x.v) - V. (pv) + pg.v 
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1 7 ~~ * 

Si les forces volumiques a distance sont uniquement celle de la gravite, on a : g = -V(gz) 
D'ou: pg.v = -(Vgz).(pv) 

Appliquons une formule du type f'g=(fg)'-fg' : (Vgz).(pv) = V.(gz.pv) - (gz)V.(pv) 



f)p 

Multiplions la conservation de la masse par (g.z) : (gz) — + (gz) V.(pv) = 

at 

dp dpgz 

D'ou (comme g est independant du temps) : pg.v = -V.(gz.pv) - gz— = -V.(pgz. v) - — - — 

at at 

Finalement, l'equation de conservation de l'energie totale (interne+cinetique+potentielle) s'ecrit: 



^ U+ a t V2 + g ^ + ^(p(u + ^v 2 + gz)v) =-V.(j q )+ r q - V.(t.v) - V.(p.v) 



Remarque. 

On peut faire apparaitre l'enthalpie: H = U + pV = U + p / p ou V est le volume massique 



gft+pg + V.(p(H + 1^) = -v4) + r, - V.(f.v) 



3.1.2) Energie interne 

Au paragraphe 2.10 nous avions derive l'equation de conservation de l'energie cinetique 
(macroscopique). 

- + V.(^pv 2 v) = -V.(pv) + p(V.v) -V.(f.v) + f:(Vv) + pg.v 
Si nous soustrayons l'equation de conservation de l'energie cinetique macroscopique (iv 2 ) a celle 



2 

resultant du premier principe de la thermodynamqiue pour l'energie totale (U+jV 2 ), nous 
obtenons l'equation de conservation de l'energie interne: 



apu 
at 


V. (pUv) = 


- v.(U 


+ r q 


- x:(Vv) 


- p(v.v) 


accumulation 


bilan des flux 
convectifs 


bilan des flux 
diffusifs de 
chaleur 


source 
volumique 
de chaleur 


energie 
mecanique 
dissipee en chaleur 
(irreversible) 


energie mecanique 
convertie de fa§on 
reversible en 
energie interne 



Remarque : le dernier terme est nul pour un fluide incompressible. 



17 II suffit en fait que ces forces derivent d'un potentiel constant dans le temps 4><->gz 
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3.2) Equation de conservation sur une particule de fluide suivie dans son mouvement . 

(derivee particulaire) 

De la meme facon que pour les equations de conservation de la masse d'un constituant et de la 
quantite de mouvement, les equations precedentes peuvent etre reecrites en utilisant la derivee 
particulaire: 

d(u + ^) _ h . — / = \ - , > 

P — = - V.(j q J + r q - V. (x.v) - V. (p.v) + pg.v 

D(u + ^ + gz) . v _ ,_ . _ , . 

P^ = " Hi) +r , - ^(Vv)-p(v.v) 

3.3) Equation d'etat 

Dans le cas general, l'energie interne massique U est une fonction de la temperature et de la 
pression (eventuellement de la composition) 

U = f{p,T} 

Pour les gaz parfaits on a: 

U = C V (T-T ref ) + U ref et H = C p (T-T ref ) + H ref 

Pour les liquides et les solides la pression a peu d'influence. En l'absence de changement de phase 
et si la capacite thermique massique : C , est constante on a: 

U = C(T-T ref ) + U ref et C p =C v =C 

Dans le cas d'un melange diphasique liquide-vapeur d'un corps pur, si x est la fraction massique de 
la phase vapeur on a: 

p = Psat {T} et U = x.U vapsat {T} + (l-x).U liqsat {T} 

Dans le cas d'un melange diphasique solide-liquide d'un corps pur, si x est la fraction massique de 
la phase liquide on a: 

T = T fu S ion et U = x.U liq {T fusion } + (l-x).U sol {T fusion } 

Si Ton a affaire a un melange dont la composition varie, l'energie interne massique est fonction en 
plus de la fraction massique des differents constituants. 

U = f{p,T,(O i } 
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A une pression et une temperature donnees, l'energie interne d'un melange n'est pas a priori egale a 
la somme des energies internes de chacun des constituants purs a ces memes pression et 
temperature, (inversement si on melange, de facon adiabatique, un litre d'eau a 20°C et un litre 
d'alcool a 20°C, le melange ne sera pas exactement a 20°C). 

On a cependant souvent affaire a des milieux dimes, dans ce cas on peut considerer que la capacite 
calorifique est constante. 

U = C V (T-T ref ) + U ref 

En ce qui concerne la masse volumique, nous avons deja indique qu'elle etait fonction, a priori, de 
la composition, de la pression et de la temperature. Dans beaucoup de problemes de circulation 
d'air avec transfert thermique (chauffage d'une piece par un radiateur par exemple) on utilise 
1' approximation de Boussinesq. Cela revient a dire que le volume occupe par un kg de fluide ne 
varie pas beaucoup avec la temperature. De ce point de vue on peut considerer que p est constant. 
Par contre il est important de tenir compte du fait que la force de gravite exercee sur un m d'air , 
soit p.g, depend de sa temperature. De ce point de vue on doit considerer que p depend de la 
temperature. On utilise pour cela une approximation lineaire autour d'une temperature moyenne 
T . 

dans toutes les equations on utilise : p = est = p = p{T } 

sauf dans le terme p.g que Ton approche par : p.g = p [l - (3(T- T )]g 
pour un gaz parfait |3 = 1/T Q 

Ceci permet, d'une facon simple, de tenir compte que de l'air chaud a tendance a monter et de l'air 
froid a descendre. 



3.4) Equation de transfert (loi de Fourier) 
3.4.1) Transfert unidirectionnel 

Nous avons vu en thermodynamique que la loi de Fourier etait l'expression, pour un milieu 
continu, du fait que le flux de chaleur va du plus chaud au plus froid et qu'il est proportionnel a 
l'ecart de temperature. Si la temperature est uniquement fonction de y alors le flux diffusif (pour la 
chaleur on parle egalement de conduction) s'ecrit : 

J -~" X dy 

On peut donner une interpretation de la conduction dans les gaz similaire a celle des transferts de 
matiere et de quantite de mouvement. Dans les solides la transmission de la chaleur peut se faire 
par d'autres mecanismes notamment par l'intermediaire du 'gaz d'electrons' dans les metaux. 
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L'agitation moleculaire (on parle aussi d'agitation thermique) des particules est d'autant plus 
importante que la temperature est elevee. Autant de particules traversent le plan perpendiculaire a 
oy de haut en bas que de bas en haut (conservation de la masse) 18 . Mais les particules qui 
descendent ont une energie cinetique plus importante que celles qui montent, d'ou un flux net 
d'energie cinetique (d'agitation moleculaire). Rappelons que pour un gaz parfait monoatomique 
l'energie interne n'est autre que l'energie cinetique des particules. Ce flux d'energie cinetique 
(d'agitation moleculaire) correspond done a un flux diffusif d'energie interne. 



A temperature elevee 



temperature basse 



On peut utiliser la loi de Fourier pour estimer le flux diffusif d'energie interne ou de chaleur dans 
la plupart des cas, en phase gazeuse, liquide ou solide. La conductivity thermique est, en general, 
fonction de la temperature et eventuellement de la pression et de la composition. 

X = f{T,p,0)J 

Si la capacite thermique a pression constante est independante de la temperature, on peut encore 
formuler la loi de Fourier en terme de gradient de l'enthalpie massique (forme plus proche de la loi 
de Fick qui est 1' equivalent en transfert de matiere) 



3T 



x ac D T 



dy H pC D dy 



si C est independant de la temperature 



3H I 

soit j q = -p.a-r— ou a = — est appele diffusivite thermique 

y p 



3.4.2) Transfert pour un champ quelconque de temperature 

Comme pour les transferts de matiere, on generalise la derive de la temperature par rapport a y par 
la notion de gradient. La loi de Fourier s'ecrit done: 



dx 
3T 

dy 
3T 

^dzJ 



et si C est independant de la temperature j = -p. a VH 



Remarque : pour les transferts de chaleur il existe trois modes de transport : la convection , la 
diffusion (ou conduction) et le rayonnement. 



La vitesse plus importante des particules a haute temperature est compensee par une densite plus faible 
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3.5) Conditions aux limites et aux interfaces 



3.5.1) Conditions aux limites 

Comme pour les transferts de matiere, les conditions aux limites se classent en 3 categories 
Dirichlet : on connait la temperature (< — > on connait les fractions massiques) 
Neumann : on connait le flux de chaleur (< — > on connait les flux des constituants) 
Cauchy : on connait une relation entre la temperature et le flux de chaleur. 



a) Condition de Dirichlet : on connait la temperature 

C'est typiquement le cas a l'entree d'un systeme. On connait la temperature de l'air souffle dans 
une serre, une chambre froide... 

C'est souvent le cas s'il y a des equilibres de phase. Si, par exemple, une paroi de conductivite 
elevee est chauffee par condensation de vapeur d'eau saturante a pression constante, la temperature 
du fluide a chauffer en contact avec la paroi sera proche de la temperature de saturation 
correspondant a la pression de la vapeur. Elle sera en fait legerement inferieure du fait de la 
resistance thermique du film de condensat et de celle de la paroi. 



Tsat{p} 



~Tsat{p} 



vapeur 



film de condensat 




fluide a chauffer 



paroi conductrice 



C'est egalement le cas si une paroi est chauffee (ou refroidie) par un systeme de regulation a 
temperature constante. (plaque electrique avec palpeur). 



b) Condition de Neumann : on connait la densite de flux diffusif de chaleur 

C'est typiquement le cas pour une paroi parfaitement isolee. La densite de flux total (convectif + 
diffusif) est nulle, par ailleurs, la vitesse est nulle pres d'une paroi, il n'y a done pas de flux 
convectif. On en deduit que la densite de flux diffusif est nulle. 

On peut souvent considerer qu'a la sortie d'un dispositif la diffusion est negligeable devant la 
convection 19 . 



II faut pour cela que la vitesse en sortie soit suffisamment importante, la tuyauterie assez longue jusqu'au 
prochain dispositif, ceci compare a la diffusion thermique moleculaire ou turbulente, v.L/a »1. Ecrire cette 
condition est suffisant si en tout point de la section de sortie, la vitesse est dirigee vers l'exterieur (absence de 
recirculation passant par cette section). 
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C'est egalement le cas si on chauffe a puissance imposee (plaque electrique sans palpeur). La 
densite de flux d'energie apportee sous forme electrique (par les resistances) est imposee et 
penetre dans le fluide par conduction (la vitesse du fluide etant nulle a la paroi) 



vitesse nulle — > transfert par diffusion (conduction) 

(sous l'effet du gradient de temperature) 



densite de flux 
connue 



domaine d' etude 
des transferts de chaleur 



L'explicitation de la densite de flux diffusif de chaleur, que nous noterons j q . en trant est tout a fait 
similaire a celle de la densite de flux de matiere (cf. § 1.5). 



Pour une paroi perpendiculaire a ox 



J q .entrant 



~ jq.x =-k 



3^ 

dx 



Dans le cas d'une paroi perpendiculaire a un vecteur n (vecteur normal unitaire oriente vers 
l'exterieur du domaine) : 



interieur du domaine 




j q .en,ran, * ~ J Q • 11 = M VT) . n 



Pour une paroi parfaitement isolee (adiabatique) on aura done 



Wt=MVT).n = => (VT).n = 



1 3T 1 3T 
pour une paroi a 45 degres : —j=— + —j=— = 

V2 ox V2 oy 



c) Condition de Cauchy : on connait une relation (affine) entre la densite de flux 

diffusif de chaleur et la temperature 

C'est typiquement le cas pour une paroi de conductivite connue (n'accumulant pas de chaleur) qui 
est en contact, du cote exterieur, avec un liquide de temperature connue et fortement agite. Sur la 
surface exterieure de la paroi, la temperature est alors proche de la temperature moyenne du 
liquide. 
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L ext 



milieu exterieur 
parfaitement agite 



^"temperature 
uniforme 




interieur du domaine 



gradient de temperature 
perpendiculairement a la surface 



parol 



Le flux traversant la paroi (par conduction) s'ecrit : 
T -T 

jq emmnt = ^ paroi ~ ~ = h - ( T ext ~ T s ) h : coefficient de transfert de chaleur de la paroi 

^ paroi 



Ce flux entre dans le domaine uniquement par conduction (diffusion) car la vitesse du fluide est 
nulle a la paroi (pas de convection). D'ou : 



h(T ext -T s ) = XVT| n 



Pour une surface perpendiculaire a ox, on a par exemple: - [ X — I = h. (T ext - T| x=0 ) 

x=0 



3T 



3.5.2) Conditions d' interface 

Dans certains cas on est amene a traiter un systeme de deux milieux continus separes par une 

interface. Cette interface peut etre la surface de separation de deux phases (liquide-vapeur, ou 
fluide-solide) ou une paroi de faible inertie thermique (cette derniere restriction ne s'applique pas 
au regime permanent). 

Contact direct entre deux milieux 

Si aucune paroi ne separe deux milieux, la temperature des deux milieux est identique a l'interface 
et le flux d'energie interne est continu (ce qui sort d'un milieu rentre dans le second). Le flux de 
masse totale est evidemment continu lui aussi. 

* X — T 

interface fixe: * (pU.v+j q ) s .n = (pU.v+ j q ) .n 

* (pv) s _.n = (pv) s+ .ri 

Si l'interface se deplace, il faut se placer dans le repere de l'interface, ce qui revient a considerer les 
vitesses relatives a celle de l'interface: v- 



* T — T 



interface mobile: * (pU.(v - Vj) + j q ) .n = (pU.(v - v ; ) + j q ) .n 
* (P(v-Vj)) ._.n = (p(v-v i )) .5 
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Prenons l'exemple de la congelation d'une couche de glace au niveau d'une plaque froide a 
temperature imposee. Nous supposerons que les capacites thermiques de l'eau et de la glace sont 
constantes. On a: 




Au niveau de l'interface : T s {z = z ; } = Tj {z = z ; } = T, 



Ainsi: X. 



fusion 

Ps( V s - V i) = Pl( V l - V i) 

P s U s (v s -v i )-^ s -^ = p 1 U 1 (v 1 - v i)-^i"^ 



_j -J - J = p s (U, - U s )Vj ou (U, - U s ) est la chaleur latente de fusion 



la difference entre la chaleur qui arrive et qui part de l'interface par conduction (terme de droite) 
sert a congeler de l'eau a l'interface ( terme de gauche) 

Remarque: la temperature est continue a l'interface mais elle n'est pas continument derivable, il 
existe cependant des derivees a gauche et a droite. 



Interface constitute par une paroi 

Dans le cas d'une paroi, si la temperature du fluide en contact direct avec la paroi est connue d'un 
cote de la paroi il suffit de calculer les transferts de l'autre cote. Si, au contraire, on ne la connait a 
priori ni d'un cote ni de l'autre de la paroi, on devra traiter le systeme constitue par les deux 
milieux et ecrire une relation de transfert a travers la paroi ainsi qu'une equation de bilan au niveau 
de l'interface. On considere en general que la paroi n'accumule pas d'energie interne (regime 
permanent ou paroi de faible inertie thermique) et on neglige en general l'epaisseur de la paroi 
d'un point de vue geometrique. Prenons l'exemple d'une paroi perpendiculaire a ox. 

temperature 
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On peut ecrire: 



3x, 



x=0- 



V 3x y 



x=0+ 



= h(T x=0 _ - T x=0+ ) avec h = 

e. 



parol 



flux arrivant a gauche = flux sortant a droite = flux traversant la paroi 

Nous avons bien deux equations scalaires, elles remplacent les 2 conditions aux limites qu'il 
faudrait ecrire pour les deux milieux pris separement (une equation scalaire pour chaque condition 
aux limites) 

3.6) Combinaison des equations - Adimensionnalisation 
3.6.1) Combinaison des equations 

Nous allons maintenant combiner les differentes equations dans le cas d'un fluide 

- de capacite thermique massique: C et de composition constantes 

- en considerant 1' approximation de Boussinesq pour la dependance en temperature 
de la masse volumique: p 

- repondant a la loi de Fourier avec une conductivity independante de la temperature 

- de rheologie newtonienne avec une viscosite cinematique independante de la temperature 

- sans apport volumique de chaleur 

- avec des conditions aux limites a flux nul ou temperature imposee 

Prenons l'exemple schematise ci-dessous et raisonnons en deux dimensions. La vitesse a l'entree 
est v , la hauteur du dispositif est L. 

/ s 1 parois adiabatiques ,. . , , 

^ M limites du domaine 






a) Interieur du domaine 
conservation de la masse 
dp 



temperature imposee 
L2 



7\ 
LI 



L 



> 



at 



+ V.(pv) = p = cste (equation d'etat) => V.v = (eq.l) 



conservation de la quantite de mouvement 

3pv 



at 



+ V.(pvv + tc) = pg © p = cste © 71 = pi - (i(Vv + (Vv) ) (loi de transfert) 

3pv - 

-r— + V.(pvv) = -Vp + (iV 2 v + pg 
at 
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Si nous utilisons 1' approximation de Boussinesq en prenant To comme temperature de reference 
pg = p g - p (3(T - T )g = -Vp gz - p (3(T - T )g 

3pv - 

=> + V.(pvv) = -V(p + p gz) + uV*v - p (3(T - T )g 

<=> ^ + V.(pvv) = - V(p) + uV'v - p (3(T - T )g 

avec p = p + p gz pression motrice (eq. 2) 

conservation de l'energie interne: 

dpU - = - 

— - — h V.(pU. v + j ) = r - x:(Vv) (loi de conservation ) 
at 

p = cste et U = C(T - T ref ) + U ref (lois d' etat ) 
j q = -XVT (loi de transfert ) r q = 

^ ^ + V.(Tv-^VT) = -^(Vv + (Vv) t ) :(Vv) 
ot pC pC V ' 

^ + V.(Tv)-aV 2 T = ^(Vv + (Vv) t ) :(Vv) (eq. 3) 
b) Limites du domaine 



Sur la section d'entree 1: nous connaissons la temperature T et la vitesse v G 



sur SI (x = 0, < z < LI) : v = 



T = T (eq.4) 



Sur la section de sortie 2: nous negligeons la conduction et nous connaissons la pression motrice 
fV 

dT 

sur S2 (x = L2, < z < LI) : — = 0, p = p (eq.5) 

Les parois 3, 4 et 5 sont adiabatiques (parfaitement isolees), la vitesse et le flux de chaleur sont 
nuls: 

- 3T 

sur S3 et S4 (x = ou L2, LI < z < L) : v = 0, — = (eq.6) 

dx 

- 3T 

sur S5 (0 < x < L2, z = L): v = 0, — = (eq. 7) 

oz 

Au niveau de la surface chauffante 6, la temperature est imposee et la vitesse est nulle 
sur S6 (0 < x < L2, z = 0): v = 6, T = T p (eq.8) 



76 



c) conditions initiales 



Nous considerons qu'initialement la temperature dans le fluide est uniforme dans le domaine et 
egale a celle du fluide entrant: T Q et que la vitesse est nulle. 

pour t=0 (quels que soient x et z): v = , T= T G . (eq.9) 



Le probleme est ainsi entierement defini, nous avons combine 

- les equations de conservation (masse, quantite de mouvement, energie interne) 

- les equations d'etat ( p = cste et U = C(T - T ref ) + U ref ) 

- les lois de transfert (fluide newtonien, loi de Fourier) 

- les conditions aux limites 

- et les conditions initiales dans un cas particulier. 



5.6.2) Adimensionnalisation 

Definissons maintenant des variables sans dimension: 

T-T n 



la temperature reduite: T* = 



l 



T -T 

p o 



- elle vaut pour le fluide entrant ou initialement present dans le domaine 

- elle vaut 1 au niveau de la paroi 

x 

- les coordonnees adimensionnelles: x* = — 

v 

- la vitesse adimensionnelle: v* = — 

v o 

En effectuant les changements de variables correspondant a cette adimensionnalisation on obtient: 
a) Interieur du domaine 
Equation de continuity 

V*.v* = (eq. 1' : conservation de la masse, eq. 1 multipliee par L/vq) 

( d d d s - 

ou V * est le vecteur 



dx * dy* dz *j 
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Equation de l'energie interne 



+ — — V*.(T*v*) - V 2 *(T*) = — lv*v*+(V*v*) ) :(v*v*j 

a A,(T -T ) 



at 



L 2 



(conservation de l'energie interne, eq.2 multipliee par L 2 /cc) 



On voit apparaitre un temps adimensionnel appele nombre de Fourier et note Fo. 

a.t t 

t* = Fo = 



L 2 L 2 / a ^~ temps caractensttque 

Considerons un milieu semi infini de vitesse nulle (les transferts se font uniquement par 
conduction), la temperature initiale est uniforme et egale T . Si a l'instant t=0, on porte la 
temperature a la surface (x=0) a une valeur Ti, alors a la profondeur L, la temperature reduite:(T- 
T )/( Ti-T ) atteint une valeur de 0,48=1/2 au bout d'une duree egale a L 2 /a 

L 2 /a est done le temps caracteristique de diffusion (temps necessaire pour que "l'information" arrive 
jusqu'a L). 

t=0 t=L 2 /a 



l l 



l0_ 




infini 




infini 



II apparait egalement le rapport (voL/a) que Ton appelle nombre de Peclet , note Pe. On peut 
l'interpreter comme le rapport entre une difference de flux convectif (vitesse v D , temperature T p ou 
T Q ) et un flux conductif (gradient de temperature (Tp-T )/L ) 



Pe = 



v n .L p.v .C.(T -T ) flux convectif 



a 



(T p -T ) 



flux diffusif 



Enfin, on voit apparaitre le groupement: Br= 



A-(T p -T ) 



appele nombre de Brinkmann. 



C'est l'ordre de grandeur du rapport entre l'energie mecanique dissipee (de fa§on irreversible) en 
chaleur et l'energie transmise a la paroi. 

* (I. v / L est l'ordre de grandeur des contraintes visqueuses, la puissance de ces forces est de 
l'ordre de (iv 2 , / L (multiplie par la surface ou elles s'exercent) 

* X(T p - T ) / L est la densite de flux de chaleur en regime permanent a travers une epaisseur 
de fluide immobile L, si la difference de temperature est T p - T (il faut multiplier ce terme 

par la surface pour obtenir le flux ). 
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L'equation de conservation de l'energie interne peut done se reecrire: 
dT* 



3Fo 



+ Pe. V* .(T * v *) - V 2 * (T *) = Br(v * v * +( V * v*) ' ) : (v * v *) (eq. 2') 



Equation du mouvement: 

dv* 



\ L J 



+ V.(v*v*) = -V* 



P-Po 
V pVo J 



+ 77tT V2 *(v :|: )- 



pcr p -T )Lg 



V V J 



On voit apparaitre un temps adimensionnel different de celui relatif a l'equation de l'energie 
interne, mais on peut remarquer que: 



= Fo.Pe 



On voit apparaitre les nombres d'Euler: 



P-Po 

V PVo J 



et de Reynolds 



v L 



V V 



g.P(T -T )L 

Enfin, on voit apparaitre le groupement: Ri = appele nombre de Richardson 

II compare les forces de poussee d'Archimede aux forces d'inertie du fluide, en effet: 

^P Vq est la pression dynamique du fluide pour une vitesse egale a v Q 
(ou encore l'energie cinetique par unite de volume) 

p g.P(T p -T )L represente la poussee d'Archimede qui s'exercerait (par unite de 

surface) sur une couche de fluide d'epaisseur L portee a la 
temperature T p et immerge dans le meme fluide a une 
temperature T D 



F=A.Lgp(3(Tp-To) 



'0 



surface:A 




L'equation de conservation de la quantite de mouvement peut done se reecrire: 



1 3v* - - 1 

— - — + V.(v*v*) = -V*Eu + — V 2 *(v*)-Ri.e T* 
Pe3Fo Re 



ou e z est le vecteur unitaire oriente selon oz (vers le haut) 



(eq 2') 
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b) Limites du domaine : 



surSl (x* = 0, 0< z*< Ll/L) : T* = 0et v* : 



surS2 (x* = L2/L, 0< z*< Ll/L) : Eu = 0, 



dx * 



= 



- 3T* 

sur S3 et S4 (x* = ou L2 / L, LI/ L < z* < 1) : v* = 0, ^ = 



3x* 



sur S5 (0 < x* < L2 / L, z* = 1): v* = , 



3T* 

dz * 



sur S6 (0 < x* < L2 / L, z* = 0): v* = , T* = 1 



(eq.4') 
(eq. 5') 

(eq. 6' et 7') 
(eq. 8') 



c) Conditions initiales 

pour t* = Fo = 0, quels que soient x * et z * : v* = 0, T* = (eq. 9') 

Les equations 1' a 8' decrivent entierement le systeme. Elles permettent de calculer revolution du 
champ des variables reduites (T*, v*,Eu) = f{x*,z*,Fo = t *} en fonction des parametres 
adimensionnels: Pe, Re, Ri, Br et de la geometrie: (Ll/L) et (L2/L) 

Remarque: 

II est possible d'obtenir directement les parametres adimensionnels (T*,x*,z*,Fo,Re,Pe,Br) dont 
depend le probleme par une analyse de tous les parametres qui entrent en jeu 
(T,T p ,T ,x,z,t,L,v ,p,a,C,(J,,g) et des dimensions fondamentales que font intervenir ces 

parametres. 

II faut avant tout remarquer que la temperature initiale et d'entree du fluide T Q peut etre prise pour 
reference. Pour une difference de temperature (T p - T G ) donnee, un decalage de T Q entraine un 
decalage identique de toutes les temperatures. 

On forme alors la matrice des dimensions formee par les variables. 





T-T 


Tp-T 


X 


z 


t 


L 


v 


p 


c 




a 


g 


m 








1 


1 





1 


1 


-3 


-1 


-1 


2 


1 


s 














1 





-1 





-2 


-1 


-1 


-2 


K 


1 


1 




















-1 











Kg 























1 


1 


1 









Attention ne pas oublier l'acceleration de la pesanteur, elle est differentes pour 
les martiens pourtant la physique est la meme pour eux. 

Nous avons 12 variables dimensionnelles en relation et le rang de la matrice des dimensions est 4, 
on peut done exprimer la relation en faisant intervenir (12-4)=8 nombres sans dimension. 
T*=f{x*,z*,Fo,Re,Pe,Br,Ri} 
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Dans la plupart des cas, la dissipation d'energie mecanique en chaleur est negligeable devant le flux 
transfere a la paroi chauffante, le nombre de Brinkmann est alors faible et n'intervient pas. Si de plus, 
les forces de poussee d'Archimede sont faibles devant les forces d'inertie du fluide, le nombre de 
Richardson est faible et n'intervient pas. On dit que Ton est en convection forcee. On a alors : 

T*=f{x*,z*,Fo,Re,Pe} 

Dans le cas particulier ou v est nul (pas de ventilateur ou de pompe), ce sont uniquement les 
forces de poussee d'Archimede qui generent le mouvement, on dit que Ton est en convection 
naturelle. Les nombres de Richardson et de Reynolds n'ont alors plus de sens pris separement. 
Par contre en prenant Ri.Re 2 on elimine vo et Ton obtient un nombre qui caracterise 1' apparition de 
circulations de fluides dues a des differences de densite. 

Tl . -+A K A n U f n g-P(V T o) L3 

II s agit du nombre de Grashof : Gr = 



3.7) Notion de coefficients de transfert entre une paroi et un fluide 

Nous avons defini precedemment un coefficient de transfert de chaleur a trovers une paroi comme 
la densite de flux de chaleur divisee par la difference de temperature de part et d' autre de la 
paroi 20 . On peut definir de la meme maniere des coefficients de transfert, notes h, entre une paroi 
et un fluide comme un flux de chaleur ramenee a une surface et a une difference de temperature 
paroi/fluide. 

On adimensionnalise ces coefficients de transfert sous la forme d'un nombre de Nusselt en les 
multipliant par une dimension caracteristique et en les divisant par la conductivite thermique du 
fluide. 

Q/S h.D 

h = Nu = D: dimension caracteristique 

T — T A, 

paroi fluide fluide 

Ecoulement autour d'un objet - temperature de surface uniforme 

Un fluide s'ecoule autour d'un objet dont la temperature de surface Ts , est uniforme. Loin de 
l'objet, la vitesse et la temperature du fluide sont uniformes (on les note Voo et Toe). On definit en 
general le coefficient de transfert sous la forme : 

Q/S Q flux total objet -> fluide 

h = 

T s - T^ S surface totale de l'objet 

On remarquera que la densite de flux de chaleur n'est pas uniforme sur l'objet (plus elevee a 
l'avant qu'a l'arriere) 



Voc 

T 

J- CX 



Nous avions egalement defini des coefficients de transfert de matiere a travers une membrane comme la densite de 
flux du constituant divisee par la difference (de part et d' autre de la membrane) de la fraction massique , de la 
concentration, de la pression partielle ... (cf § 1.5) 
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On peut ecrire l'ensemble des equations regissant ce probleme et (en principe) trouver le champ 
des temperatures en fonction de la vitesse loin de l'objet, de la viscosite, des temperatures et T s> 
etc. 

Le flux de chaleur qui sort de l'objet a travers sa surface S, entre dans le fluide par conduction 
(vitesse du fluide nulle a la surface), il s'ecrit done : 



Q = \ -XVT.ndS (X est celui du fluide) 




Connaissant la temperature en tout point, on peut en deduire le gradient pres de la surface et done 
le flux de chaleur sortant. Finalement on peut theoriquement exprimer le coefficient de transfert en 
fonction de la vitesse loin de l'objet, de la viscosite, etc. 

Sous forme adimensionnelle, avec les hypotheses de convection forcee en regime permanent, on 
trouve : T*=f{x*, Re, Pe} d'ou l'on tire Nu=f{Re,Pe} 21 . 

On peut bien sur remplacer, dans cette relation, le couple (Re,Pe) par (Re, Pe/Re). Le rapport 
Pe/Re est appele nombre de Prandtl, il a l'avantage de ne dependre que des caracteristiques du 
fluide et non pas de l'ecoulement. II represente le rapport entre le coefficient de diffusion de la 
quantite de mouvement (par les forces visqueuses) et le coefficient de diffusion de la chaleur. 

Pr = 

Ainsi en convection forcee on peut ecrire : Nu = f { Re, Pr }. Les relations de ce type peuvent etre 
obtenues par une resolution analytique des equations dans des cas tres simples, par resolution 
numerique pour des geometries complexes en ecoulement laminaire, mais le plus souvent elles 
sont identifiers par des experiences. On parle de correlation entre nombres sans dimension. 22 

Theoriquement la correlation obtenue n'est valable que pour une geometrie donnee, par exemple 
une sphere isolee, et ceci pour une temperature de surface uniforme, mais l'ingenieur l'utilisera 
neanmoins pour calculer le refroidissement de pommes (^spheres) dont la temperature de surface 
n'est pas egale en tout point lors du refroidissement. II faut etre conscient que le resultat n'est alors 
qu'une approximation. 

Ecoulement autour d'un objet - densite de flux uniforme 

Supposons maintenant que la densite de flux de chaleur est imposee de facon uniforme sur la 
surface, dans ce cas la temperature de surface n'est pas uniforme. On prend alors, dans la 
definition du coefficient de transfert, la temperature de surface moyenne. 



=1 



T dS/ S h = = — — q : densite de flux imposee a la surface 



L'expression de la correlation Nu{Re,Pr} sera a priori differente que l'on soit a temperature de 
surface uniforme ou a densite deflux de chaleur uniforme. 



21 La position x* n'apparait pas puisqu'il s'agit d'un coefficient de transfert moyen sur toute la surface de l'objet. 

22 Les expressions C x {Re }sont des correlations similaires en mecanique des fluides 



82 



Ecoulement dans une conduite - flux de chaleur uniforme 

Un fluide s'ecoule dans un tube raccorde a un grand recipient. Admettons que Ton soit en regime 
laminaire 23 , le profil parabolique de vitesse ne s'installe pas immediatement mais seulement apres 
une certaine longueur dite d'etablissement. Prenons le cas ou, sur cette zone d'etablissement 
hydrodynamique la paroi du tube est parfaitement isolee, alors qu'au dela une densite de flux de 
chaleur uniforme q est imposee a la paroi. 



Profil de vitesse . 


. Profil de temperature 

llllllllll 
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Lorsque Ton avance dans la zone de chauffage, pres de la paroi, une couche de plus en plus epaisse 
de fluide va etre chauffee. On parle de couche limite thermique, une theorie tres utilisee mais qui 
sort du cadre de ce cours permet de prevoir l'epaisseur de cette couche. Lorsque la couche limite 
atteint le centre du tube, on dit que le regime thermique est etabli. Le profil de la temperature reduite 

T* = (T{r,z} -T{r = 0,z}) / (T p {z} -T{r = 0,z})ne depend alors plus de z. 



On peut definir un coefficient de transfert local a partir de la densite de flux de chaleur a la paroi, de 
la temperature de la paroi a la position axiale z consideree et d'une temperature moyenne du fluide 
sur la section : le plus souvent on prend la temperature de melange T{ z} . 

a [(T.vJdS f R T{r,z}v{r}.27ir.dr 

l0CaltZj ~T{z}-T{z} aV6C f vdS tubec^ai. V 

J s z 

On peut ecrire l'ensemble des equations regissant ce probleme et calculer (en principe) le champ de 
temperature, notamment T p {z} etT{z}, en fonction des donnees du probleme : debit, densite de 

flux de chaleur, etc. On peut done en deduire le coefficient de transfert. Sous forme 
adimensionnelle, avec les hypotheses de convection forcee en regime permanent, on trouve : 

Nu local = fi { z/D, Re, Pr } fi coefficient local & densite de flux uniforme 

L'interet d'utiliser la temperature de melange est que, dans certaines conditions, pour l'equation de 
bilan d'energie interne, tout se passe comme si la temperature du fluide etait homogene, egale a T . 
En effet, en regime permanent , pour un liquide de capacite thermique volumique pC independant 

24 

de la temperature et si la diffusion axiale est negligeable , le bilan sur une longueur elementaire dz 
de la conduite s'ecrit: 

dff pCT.v z .27tr.dr) = q.27tR.dz => pVC.dT{z} = h{z}.(T p {z} - T{z}).27lR.dz 



On supposera egalement pour la suite que la viscosite du fluide est independante de la temperature ce qui est 

loin d'etre le cas pour des produits alimentaires. 

On neglige egalement la dissipation d'energie mecanique en chaleur. 
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On peut egalement definir un coefficient de transfert moyen sur l'ensemble d'un tube. 
h moye „{U = | o L h local {z}dz/L Nu moyen = f 2 { L/D, Re, Pr } f 2 ^ 

Les expressions des correlations sont bien sur differentes selon que Ton est en regime laminaire ou 
turbulent et qu'il s'agit d'un Nusselt local ou moyen. 

Ecoulement dans une conduite - temperature de paroi uniforme 

On peut faire un raisonnement similaire. Les correlations pour les coefficients de transfert (locaux 
et moyens) seront, a priori, differentes que Ton soit a temperature de surface uniforme ou a 
densite de flux de chaleur uniforme. 

Nu local = gi { z/D, Re, Pr } gi ( coefficient local & temperature uniforme ) 

Nu 

moyen — 

g 2 {L/D, Re, Pr } g 2 * gl g 2 *f 2 

De nombreuses correlations du type: Nu=f{z/L,Re,Pr} existent dans la litterature pour differentes 
geometries. Verifiez toujours s'il s'agit de coefficients locaux ou moyens. D'une facon generate il 
faut toujours se reporter a la relation correspondante entre flux (a definir precisement) et ecart de 
temperature (a definir precisement). Par ailleurs verifiez toujours les conditions de validite des 
correlations (gamme de Reynolds et de Prandtl). 

Pour des problemes concrets, l'ingenieur est amene a utiliser des correlations dont les hypotheses 
ne sont pas exactement verifiees. II prendra cependant soin de se rapporter a la situation la plus 
'approchante' de la litterature et gardera en memoire que le resultat est lui meme approche. 

Tout l'art de l'ingenieur consiste a faire des approximations raisonnables a defaut d'etre 
toujours raisonnees, cet art ne s'acquiert malheureusement que par la pratique. 

Transferts conjugues fluide/solide 

De nombreux problemes concrets portent sur le refroidissement ou le rechauffement d'un solide 
par l'intermediaire d'un fluide. On ecrira alors les equations de transfert de chaleur en regime 
transitoire dans le solide (uniquement par conduction) avec une condition aux limites faisant 
intervenir le coefficient de transfert entre la surface du solide et le fluide. 
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3.8) Equation de conservation sur un domaine macroscopique 



En utilisant la formule de Leibnitz et le theoreme de Green-Ostrogradskii, on peut integrer 
l'equation de conservation de l'energie totale dans le cas ou les forces exterieures sont celles de 
gravite 



25 



3p(u + ^v 2 + gz) 



at 



+ V. (p(u + ^v 2 + gz)v) = - V.(j q ) + r q - V. (tc.v) 



On obtient: 



dt 



J v p(u + iv2 + gz)dV + | s p(u + ;v2 + gz)(v-v s )iidS = 



accumulation d'energie bilan des flux convectifs 
totale d'energie totale 

-\M dS + U dV "j s P^dS -J s v.(^.n)dS 

bilan des flux flux de chaleur puissance des puissance 
diffusifs de apporte au sein forces de des forces 

chaleur du volume pression visqueuses 

(systeme pouvant avoir des frontieres mobiles se deplacant a la vitesse v s ) 



Si les frontieres du systeme sont fixes et en negligeant les flux diffusifs sur les sections d'entree et 
de sortie, on obtient: 



dt 



| v p(u + ^v 2 + gz)dV + j" s p(u + pV + ^v 2 + gz)v.ndS Q + W e . 



accumulation d' energie 
totale 



bilan des flux convectifs flux de chaleur apporte puissances des 
d enthalpie totale au niveau des parois et forces exterieures 

au sein du volume sauf pression et gravite 



Si on a affaire a d'autres forces derivant d'un potentiel 4> , remplacer gz par 4> 
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Pour un systeme de frontiere fixe, ayant une section d'entree et une section de sortie et en 
negligeant les flux diffusifs (de chaleur et de quantite de mouvement) sur les sections 
d'entree/sortie, on peut ecrire: 




E tot est l'energie totale (interne, cinetique et potentielle) du fluide contenu dans le domaine : 

E tot ={ v p(U + |v 2 + gz)dV 

H tot est l'enthalpie massique totale : H tot = U + pV + 1 v 2 + gz 

V est le volume massique (V = 1/ p) 

(en general p et V sont considered comme uniformes sur les sections d'entree/sortie) 

i v 2 est l'energie cinetique de melange : v 2 = £ pv 2 ((v - v s ) .n)dS^m 

- si la direction de la vitesse est identique sur toute la section d'entree et que 

p est uniforme sur cette section, alors v 2 = j" s v 3 dS ^/ j^v-dS 

- si la vitesse est uniforme on a simplement : v 2 = v 2 

Q est le flux total de chaleur apporte au fluide contenu dans le systeme par les parois 

et au sein du volume (le flux de chaleur convectif sur les section d'entree/sortie est 
pris en compte par ailleurs) 

W ext est la puissance des forces exterieures (sauf celles dues a la gravite et a la pression 
sur les section d'entree et de sortie), il s'agit typiquement de la puissance mecanique 
apportee par les pales d'une pompe ou d'un ventilateur (ou de la puissance mecanique 
sortant du systeme dans le cas d'une turbine : terme negatif dans ce cas) 

Remarque: en regime permanent et si Ton peut negliger les variations d'energie cinetique et 

potentielle entre l'entree et la sortie ,ona: AH = (q + W ext ) / rh 

s— e 

Rappel: une expression plus generale (frontiere mobile) a ete fournie en thermodynamique (l er principe) 



86 



4) Introduction a l'approche filaire (modele ID) 



Beaucoup de situations peuvent etre modelisees de facon simplified par une approche 
unidirectionnelle bien que la geometrie reelle du systeme soit tridimensionnelle. Exemple : un fil 
est bien tridimensionnel mais pour des calculs mecaniques il suffit en general de raisonner sur 
revolution de la tension (integrate des contraintes de traction sur la section transversale) le long de 
la fibre neutre (lieu des barycentres des sections transversales) du fil. 

L'exemple que nous traiterons est celui d'un echangeur de chaleur tubulaire en regime turbulent. 
Nous aborderons d'abord la modelisation des transferts de chaleur. Bien d'autres cas relevent de 
l'approche filaire : ecoulements dans les canaux et rivieres, ecoulements et transferts dans les 
vaisseaux sanguin, ecoulements et transferts dans l'anse de Henle d'un nephron ... et toutes les 
equations de conservation (masse, masse d'un constituant, quantite de mouvement, ...) peuvent 
etre formulees de facon analogue. 

On considere 1' echangeur represente ci-dessous 



vapeur saturee a lbar 

I 




Considerons le bilan d'energie totale sur un troncon de tuyau compris entre x-Ax/2 et x+Ax/2. Si 
Ton neglige la diffusion axiale, on a (bilan macroscopique sur un systeme a une entree et une 
sortie) : 

^+ Am(H + |v 2 + gz)=Q + W ext 
dt ~ v 2 ' 



Dans le cas du chauffage d'un liquide, on peut negliger les variations d'energie de pression, 
cinetique et potentielle devant les variations d'energie interne : 

AU = CAT ordre de grandeur des variations sur lm : 4000 J.kg' 1 .K 1 * 10 K = 40000 J.kg" 1 
ApV ordre de grandeur des variations sur lm : 10 4 Pa. 10" 3 m 3 kg -1 = 10 J.kg" 1 

-v 2 ordre de grandeur :-\Q 2 (m's" 1 ) 2 = 50 J.kg 1 
2 2 

Agz ordre de grandeur des variations sur lm : 9.8 lm's" 2 . lm ~ 10 J.kg" 1 
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Par ailleurs, il n'y a pas d'apport d'energie mecanique (autre que par les forces de pression sur les 
section d' entree- sortie et de gravite) et si le fluide est de masse volumique constante, la masse 
contenu dans le systeme (troncon de tube rigide) est constante et 1' equation se resume a : 

- dTixi 
mC^^+mCAT = Q 

dt s - e 

Dans cette equation, en toute rigueur, la temperature dans le premier terme est la temperature 
moyenne surfacique sur la section transversale a la position x alors que dans le deuxieme terme, il 
s'agit de la temperature moyenne de melange . Dans le cas turbulent qui nous interns se ici, la 
temperature, comme la vitesse, est presque homogene sur une section transversale (a cause de 
l'importante diffusion turbulente radiale) sauf dans la couche limite a proximite immediate de la 
paroi. Aussi peut-on assimiler ici ses deux temperatures; il n'en est pas de meme en ecoulement 
laminaire. 

Exprimons maintenant la masse, le debit massique et le flux de chaleur a la paroi: 
m=pSAx 

rh = pSv ouv = | v x dS /S = V/S est la vitesse debitante 

Q=q)q %Ax ou (p q est la densite de flux de chaleur a la paroi (W.nr 2 ) et % le perimetre 
pour un tube circulaire (entierement rempli) %=7tD 

Divisions l'equation par Ax et passons a la limite lorsque Ax tend vers : 

- dT - dT 

pSC — + P SvC — = <p qX 
dt dx 

Remarque dT{x}/dt a ete reecrit dT/dt pour rendre compte que T depend de deux variables x et t et 
qu'il s'agit ici de la variation de T par rapport a t a x constant (les autres variables d'espace 
n'interviennent pas car T est en fait une moyenne sur la section). 

La densite de flux de chaleur a la paroi ne peut pas etre calculee dans cette approche ID par la loi 
de Fourier (qui ferait intervenir la variation radiale de la temperature au niveau de la paroi), elle 
sera exprimee a partir d'un coefficient de transfer! de chaleur entre la paroi et le fluide (voir 
paragraphe 3.7) qui lui-meme sera calcule d'apres une correlation tiree de la litterature pour le cas 
considere (ou s'en rapprochant le plus). 

(p q =h(T P -T) ou h est le coefficient de transfert de chaleur entre la paroi et le fluide 

Le coefficient de transfert depend des proprietes et de la vitesse du fluide. Sous forme 
adimensionnelle, en convection forcee, cette relation relie un nombre de Nusselt aux nombres de 
Reynolds et de Prandtl. 

Nu=f{Re,Pr}o JfUfkl) 
L 1 X L v oc J 
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Pour un ecoulement turbulent dans un tube, on peut citer notamment la correlation proposee par 
Dittus et Boelter : 

Nu=0.024 Re 8 Pr 4 en situation de chauffage pour 0.7<Pr<120, 2500<Re<1.2 10 5 et L/D>60 

(d'autres correlations existent pour le refroidissement, en ecoulement laminaire, en regime non etabli...) 

En regime permanent, en utilisant la notion de coefficient de transfert, on obtient : 
pSvc||- = xcp q =Xh(T p -T) 

En regime permanent pour un fluide incompressible, on peut par ailleurs ecrire le bilan de quantite 
de mouvement ou d'energie mecanique sur un troncon de tube (bilans equivalents dans ce cas): 

c d P* c 1 , A P* i 4S 

S-f- = -%x =-%Fa-pv 2 o - — ^r- = ^ f = 4 Fa avec D h = — 

dx 2 Ipv 2 — 1 

2 D h 

On voit que le nombre de Fanning (Fa) ou le coefficient de perte de charge lineaire (kf) 
correspondent a une sorte de coefficient de transfert qui, rappelons le, dependent du nombre de 
Reynolds. 

Nous mettrons en oeuvre cette approche pour dimensionner un echangeur de chaleur et la pompe 
qui l'alimente au cours d'un TD. 

L' approche filaire est utilisee dans de nombreux domaines (ecoulement dans les canaux et riviere, 
coup de belier..); sont indiques ci-dessous les equations de transport de la masse, de la masse d'un 
constituant et de la quantite de mouvement pour des situations assez generates (masse volumique 
et section variable) avec comme hypothese principale que les variables sont uniformes sur une 
section transversale (sauf a proximite des parois). 

3pvS 3pvSv ,3p 
dt ' dx ~ XpX b dx~ 

Les deux premiers termes des equations sont 1' accumulation et le bilan des flux convectifs. 

Pour la masse totale le membre de droite est un apport de masse au niveau des parois (si elles sont 

poreuses par exemple) 

Pour la masse d'un constituant le membre de droite comporte un terme de diffusion effective 
axiale (eventuellement du a la turbulence), un terme d' apport au niveau des parois et un terme 
source volumique (du a des reactions chimiques) 

Dans le bilan de quantite de mouvement, la courbure du systeme (canalisation ou canal) est 
negligee. Les termes sources sont les frottements en paroi ainsi que le bilan des forces de pression 
et de gravite qui sont regroupees dans la pression motrice. II faut considerer le perimetre mouille 
pour un canal partiellement rempli. 
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RECAPITULATIF Que nuln'entre ici s'iln'est geometre Platon (grave au fronton de I'Academie) 

Equations de conservation sur un element de volume fixe • (ce point marque ce qui est essentiel) 



Accumulation 

dp 
dt 


Bilan des flux convectifs et diffusifs 
(sortie - entree) 

+ V.(pv) = 


Terme source 




3pco ; 

at 


+ 


V. (pcOiV + l) 




r ; 






dpv 

at 


+ 


V. (pvv + ^) 




P-g 


- Vp 




3pU 

at 


+ 


V.(pUv + j q ) 






- p(v.v) - f : 


(Vv) 


dt 


+ 


V.(pSv + j q /T) 






+ o s avec o s 


>0 



Equations de conservation sur un element de masse fixee suivi dans son mouvement 



Accumulati on Bilan des flux diffusifs 

(sortie - entree) 



D(l/p) 
Dt 

Dco ; 



Terme source 
(V.v) 



Dt 



Dv 



DU 



Dt 

PS 
Dt 



v.j, 

V.( j q /T) 



P-g " Vp 



r q - p(V.v) - x:(Vv) 



r q /T + a s avec a s > 
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Equations d'etat • 



p = f{p,T,co i } 
U=f{p,T,to i } 

Equations de transfert (exemples) • 
Ja = -P d a B Vco a loi de Fick 

x = -pv( Vv + (Vv) ' ) fluide newtonien incompressible 
j q = -XVT (= -p.aVH si C p = cste ) loi de Fourier 



Equations aux limites 















= connu, 




-pDjVWj.ii = connu 


ou fjo);, VcQi.n} = 


V 











I 


connu, 




-XVT.ii = connu, 


ou f{T, VT.ii} = 
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Equations de conservation sur un domaine macroscopique 



cas general 

accumulation bilan des flux convectifs bilan des flux diffusifs terme 

a travers la limite S a travers la limite S source 





du domaine 


du domaine 


d f 

— pdV 
dt^ 


+ J s Pl v - v s)-ndS 





d f 

— p©=dV 


+ j^fv-vJ.ndS 


= - f 1.5 dS + f r.dV 

Js Jl Jv 1 


d f 

— I pvdV 
dt J v^ 


+ J s p(v(v-v s ))iidS + J s p n.dS 


= -J^-ads + J v pgdv 


^ J v p(x a v)dV + J s p(((x a v)( v - v s )) ii dS + J s p(x a n) .dS = - J s x a (X. n) dS + J v p(x a g)dV 


( reg.permanant" 
Vp = cons tan t , 


Jp(^v 2 + gz)(v- v s ).iidS + j" pv.ii.dS 

s 


= -j(-x:Vv)dV+ W ext 

V 




+ j^E tot (v-v s ).ndS + J s pv.iidS 


= "U- HdS + J/ q dvJv.(f.n)dS 26 


d r ~ 


+ j s pS(v-v s )iidS 


= "I Y- n dS + [,^dV + {a s dS o s >0 




— 


Relation obtenues par : 




di>\ 


Ostrogradskii : V cp{x,y,z,t} e R 
^{cp.dV = |^.dV + |(p.v s .ndS 




^ ^^- S 








Leibnitz : V (j){x, v,z,t} e R 3 
|(V.cp)dV = J(9.n)dS 

V s 


11 tot _ U + 2 


v 2 + gz 
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Equations de conservation sur un domaine macroscopique 

cas d'un systeme ayant une section d'entree et une section de sortie 

en negligeant les flux diffusifs sur ces sections 

accumulation bilan des flux convectifs bilan des flux terme 

sur les section d'entree a travers les parois source 
et de sortie 



dm 
~dt~ 



+ Am 

s-e 



dm ; 
~dt~ 



+ A(mu) i ) 



m ; 



+ 



dP 

dtT 

dfto 
dt 



+ A(rhv + pS) 



- rsaufpressione-, 
ext V et gravite 



+ mg 



A(m(xAv) + p(xAS)) = M . ext ( sauf — ">j • n„.iN, 



'reg. permanent^ 
v p = CSte ) 



Av(p + i P v2 + pgz) = Wj^SI " % (frontiere fixe) 



dE, 



dt 



A(m(H+ ' v 2 + gz) 



= Q + W e 



(frontiere fixe) 



dS 
dt 



+ 



A(mS) 



= Z TjT- + Ps avecP s >0 

k L k 

temperature de 1' interface oil intervient le flux de chaleur Qj, 




Bilan sur un composant / generalisation a plusieurs e/s 
Bilan sur un ensemble complexe 











V 




v... J 


s 
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Produit d'un scalaire et d'un vecteur 



s U = u s = 



Gradient d'un champ scalaire ( u — > V ) • 




Vco ; 



ax 



dy ) 



Produit scalaire de 2 vecteurs 



u.v = v.u = 



( U x> U y) 



= u x v x +u y v y 




composantes = evolution de C0i selon x et y 
Vcflj = vecteur 'plus grande pente' de (fy {x,y} 

Divergence d'un champ vectoriel • 

3pco i v x apco ; v 

bilan (sortie - entree) des flux a travers les facettes 
(divise par le volume de 1'element) , ici il s'agit du 
flux convectif du constituant i 



Produit tensoriel de 2 vecteurs ► Gradient d'un champ vectoriel 



u v = 



( V x . V y) 



Vv x u x v y ^ 



U„V. U y V y y 



Vv = 



dw^ 



dx dx 

3v, dv v 



1+b 



1 




1+a 



dy dy 

composantes = evolution de v x et v y selon x et y 
1> = ^(Vv + (Vv) 1 ) = 



^a 



U bj 

tenseur des vitesses de deformation 



Produit d'un vecteur et d'une matrice - Divergence d'un champ tensoriel 



^71 71 ^ 

xx xy 



U.7C = 



V^yx ^yy/ 



V.7C = 



3tt_ dn yx dn xv diz 



■ + 



xy 



yy 



( U x > U y) 



dx dy dx dy J 









~~ 

/ ^xx 



(7txy =^yx ) 



Norme d'un vecteur 

su.u = (u.u) s = (u x + u y ) s 



bilan des forces (=flux de quantite de mvt) 
a travers les facettes /volume de 1'element 

composantes = bilan des forces selon x et selon y 

■> Laplacien = Divergence du gradient • 

- 2 3 2 co i 3 2 co i r i 

V 0^= 2 + 2 = 'courbure' de C0;(x,y) 
ox dy 
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On obtient les expressions developpees des equations a partir des relations contractees en appliquant des regies de calcul 
matriciel . Pour ne pas surcharger les notations, en mecanique, on n'indique pas si les vecteurs sont transposes ou non. 
Pour savoir si un vecteur doit s'ecrire en ligne ou en colonne, on regarde si le resultat doit etre un scalaire, un vecteur ou 
une matrice ; on dit que Fordre tensoriel vaut respectivement 0, 1 ou 2. 



L' ordre tensoriel d' un produit simple est la somme des ordres tensoriels des facteurs 

Vp: 



p 



3p/3y 



(vecteur) (scalaire) — > (vecteur) 
1 1 



Vv = 



VV : 



3/3y y 
fv > ( W 



v v yy 





3v x 


*0 




3x 


3x 




3v x 


3v y 




,3y 




f 




A 






V V 


V V 

V * y 





(vecteur) (vecteur) — > (matrice) 
1 1 2 



(vecteur)(vecteur) — > (matrice) 
1 1 2 



L' ordre tensoriel d' un produit scalaire est la somme des ordres tensoriels des facteurs minore de 2 (produit contracte) 
V.v = 



(d/dx, 3/3y) 



3v v 3v v 



- + - 
3x 3y 



(vecteur). (vecteur) — > (scalaire) 
1 1 = (l + l)-2 



V.Jt: 



7t 7t 

xx xy 

7t 7t 

V V yyy 



37t xx 37t yx 37t xy 37i yy ~) (vecteur) .(matrice) — > (vecteur) 
~3x~ + ~3y~'~3x~ + ~3y~) 1 2 l = (l + 2)-2 



(3/3x, 3/3y) 

L' ordre tensoriel d' un produit doublement contracte est la somme des ordres tensoriels des facteurs minore de 4 



1 j 



AlB =trace(A.B) = ^^a ij .b jl 
Applications 

v.(p!) = 



V. (Vv) 1 



V. (Vv) = 



i) 


(o 3 


f3p 






^3x 


3 




^3v x 






3x 


3y 


3y 


• 


3v y 


3v y 






y 3x 


"37, 






3v x 


3v/ 






3x 


3x 




» 


3v x 


3v y 




V 


3y 





(matrice) '. (matrice) — > (scalaire) 
2 



Vp 



2 


= 


]-(£• 


dp] 




3yJ 



3x 2 3y3x 3x3y 3y 2 



\ 


r d 


r dw JL 




3 


^ v x 




/ 


v 3x 


K 3x 


3y J' 


3y 


^ 3x 





: [V(V.V 



3 2 v Y 3 2 v 
+ 



3x 2 3y 2 3x 2 3y 2 



3 2 v 3 2 v i / 



Pour des vecteurs et des tenseurs, la derivation d'un produit : (fg)'=f g+fg' s'ecrit: 
(on pourra verifier facilement l'homogeneite de ces relations: scalaire ou vecteur) 

V(sv) = v(Vs) + s(Vv) 
V.(sv) = (Vs).v+s(V.v) 
V.(vw) = w(V.v) + (v.Vw) 
(v. V)v = X - V(v. v) + (V A v) A V 



V.(c.v) = v.(v.aJ + a:(Vv) 
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Principaux nombres sans dimension (on note D la longueur caracteristique) 



flux convectif / flux diffusif 



resistance interne/resistance externe 



Peclet matiere : Pe = 



vD 



D 



Peclet thermique : Pe = 



AB 

vD 
a 



vD 

Peclet de quantite de mouvement = Reynolds : Re = 

v 



Biot matiere : Bi = 



k(D/2) 
D 



Biot thermique : Bi = 



A/solide 

h(D/2) 

bolide 



rapports de flux diffusif coefficient de transfert adim. 



temps adimesnionnel 



Schmidt : Sc 



D 



AB 



v //C„ 
Prandtl:Pr = — = — *- 
a A 

Sc a 
Lewis : Lc = — = 



Pr D 



AB 



Sherwood : Sh = 



kD 



D 



AB 



, xt hD 
Nusselt : Nu = 

X 



Fanning : Fa 



parol 



/TV 2 



Fourier matiere : Fo = ^ ' 



(D/2Y 



Fourrier thermique : Fo 



~~ t* = — ou t* = 

4 D V 



at 



(D/2) 2 



rapport de forces 



Reynolds : Re 



vD /?vD inertie 



jU visqueux 



Euler:Eu = ^P 1 ^= 
pv 2 inertie 



Richardson : Ri 



. _g/?ATD pousseed'Archimede 



v 2 inertie 

.3 



ej3 ATD 3 gJ3 ATD 3 ^ 
Grashof = Ri.Re 2 Rayleigh = Gr.Pr = ^tL 



Froude : Fr = 



V 2 
inertie 



va 



Newton : Ne = 



D g pesanteur 
F 



c _ Fx/S frontale frottement 



- 2 pv 2 



Weber: We = 



pv 2 D 2 

pv 2 D 

o tension superficielle 



inertie 



inertie 
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Nomenclature (non exhaustive) 

A surface d'echange 

C p capacite calorifique a pression constante 

C v capacite calorifique a volume constante 

d,D diametre 

D AB diffusivite 

E tot energie totale (interne, cinetique, potentielle) 

F force 

g acceleration de la pesanteur 

h hauteur 

H enthalpie massique 

j densite de flux diffusif 

(likg.m-V j q :W.m- 2 ) 

k constante de Boltzmann 

k coefficient de transfert de matiere 

L longueur caracteristique 

m masse 

rh debit massique 

M masse molaire 

M ext moment des forces exterieures 

n nombre de molecules par unite de volume 

n vecteur normal a la surface oriente vers l'exterieur 

P quantite de mouvement du fluide contenu dans le systeme 

Q flux de chaleur 



m 2 ) 

J.kg-i.K" 1 ) 

J.kg-i.K" 1 ) 
m) 

m 2 .s _i ) 
J) 

;n) 

m.s" 2 ) 
m) 

^-kg" 1 ) 

kg.m^s-^K" 1 ) 

m.s"l) 

m) 

;k g ) 

kg.s-1) 
kg.mol'l) 
N.m) 
m" 3 ) 

kg.m.s"!) 

;w) 
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R 


constante des s?az narfaits 


(J K"l mole" 1 ) 


r 


rayon 


(m) 


S 


surface (limite du systeme) 


(m 2 ) 


T 


temperature 


(K) 


u 


enersHe interne massiaue 




V 


vitesse 


(m.s" 1 ) 


V 


volume 


(m3) 


V 


debit volumique 


(m 3 . s" 1 ) 


w 


puissance 


(W) 


w d 


puissance mecanique dissipee en chaleur (W^ / V: perte de charge) 


a 


diffusivite thermique 


(m^s" 1 ) 


(3 


coefficient de dilatation volumique 


(K-l) 


e 


porosite 




$ 


potentiel 


(m 2 .s~ 2 ) 


r\ 


densite de flux total (convectif + diffusif) 


(rvkg.nrV, 


X 


conductivity thermique 


(W.m-i.K" 1 ) 




viscosite dynamique 


(Pa.s) 


V 


viscosite cinematique (|l/p) 


(m 2 .s-!) 


% 


tenseur des contraintes totales 


(Pa) 


n 


moment de quantite de mouvement du fluide contenu dans le systeme (kj 


e 


angle 




p 


masse volumique 


(kg.m- 3 ) 


T 


tenseur des contraintes de deformation 


(Pa) 


COi 


fraction massique 


(kg constituant , 



i2 c-b 
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